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5.20.  The  amplitude  and  phase  screen  are  uncorrelated.  The  mean  square 
error  in  the  tilt  estimation,  Tmaa  and  the  mean  square  tilt,  Tms,  are 
plotted  versus  LqAo  for  Dfro  =  10,  L/Lq  —  10000,  and  \/Lo  = 

5.0  X  10“®.  The  Tmse  is  calculated  twice,  once  using  the  phase  and 
amplitude  screen  case  and  once  for  the  phase  screen  only  case.  The 
Tms  is  calculated  by  computing  the  average  slope  of  the  phase  over 
the  screen.  Dp/Lg  =  1.5  and  D^/Lg  =  1.5 .  5-24 

5.21.  The  amplitude  and  phase  screen  are  correlated  by  the  cross  correla¬ 
tion.  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  L/Lg 
for  D/rg  =  10,  Lg/rg  =  200,  and  X/Lg  =  5.0  x  10“®.  Strehl  ratios  are 
plotted  for  the  compensated  screen  with  phase/amplitude  centroid, 
SRuTAPt  compensated  screen  with  actual  phase  tilt,  SRut,  compen¬ 
sated  screen  with  phase  centroid,  SRutp,  and  uncompensated  screen. 


SR.  Dp/Lg  =  D^/Lg  =  Df/Lg  =  D^/Lg  =  1.5 .  5-25 

5.22.  The  amplitude  and  phase  screen  are  uncorrelated.  Strehl  ratios 
of  the  simulated  wavefront  are  plotted  versus  L/Lg  for  D/vg  =  10, 


Lg/rg  =  200,  and  X/Lg  =  5.0  x  10~*.  Strehl  ratios  are  plotted  for  the 
compensated  screen  with  phase/amplitude  centroid,  SRutaPj  com¬ 
pensated  screen  with  actual  phase  tilt,  SRuti  compensated  screen 
with  phase  centroid,  SRutp,  and  imcompensated  screen,  SR.  Dp/Lg  = 

1.5  and  DpLg  =  1.5 .  5-25 

5.23.  The  amplitude  and  phase  screen  are  correlated  by  the  cross  corre¬ 
lation.  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus 
Lg/vg  for  D/rg  =  10,  L/Lg  =  10000,  and  X/Lg  =  5.0  x  10“®.  Strehl 
ratios  are  plotted  for  the  compensated  screen  with  phase/amplitude 
centroid,  SRutap,  compensated  screen  with  actual  phase  tilt,  SRut, 
compensated  screen  with  phase  centroid,  SRutp,  and  uncompensated 
screen,  SR.  Dp/Lg  =  Dp/Lg  =  D^/Lg  =  Dp/Lg  =  1.5 .  5-26 
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5.24.  The  amplitude  and  phase  screen  are  uncorrelated.  Strehl  ratios 
of  the  simulated  wavefront  are  plotted  versus  Lo/ro  for  D/vo  =  10, 

LjLo  =  10000,  and  A/io  =  5.0  x  10“®.  Strehl  ratios  are  plotted 
for  the  compensated  screen  with  phase/amplitude  centroid,  SRutaP} 
compensated  screen  with  actual  phase  tilt,  SRut,  compensated  screen 
with  phase  centroid,  SRutp,  and  uncompensated  screen,  SR.  Dp/Lo  = 

1.5  and  D^Lo  =  1.5 .  5-26 
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Abstract 

A  new  atmospheric  turbulence  screen  generator  is  developed  for  use  in  performance 
calculations  of  adaptive  optics  systems  valid  over  a  wide  range  of  atmospheric  turbulence 
parameters.  The  screen  generator  accounts  for  diffraction  effects  caused  by  weak  turbu¬ 
lence  and  incorporates  the  phase,  amplitude,  and  cross  statistics  of  the  perturbed  optical 
field.  The  wavefront’s  phase  and  amplitude  perturbations  are  taken  from  the  correlation 
functions  developed  by  Lee  and  Harp  and  the  cross  correlation  of  the  phase  and  ampli¬ 
tude  derived  in  this  thesis.  The  screen  generator  uses  a  modal  representation  to  perform 
a  Fourier  series  expansion  of  the  wavefront  phase  and  amplitude  over  a  square  area.  The 
phase,  amplitude,  and  cross  power  spectral  densities  determined  from  the  correlation  func¬ 
tions  are  used  in  the  Fourier  series  expansion.  The  mean  square  value  and  the  structure 
functions  of  the  phase,  amplitude,  and  cross  terms  are  calculated  to  within  1%  of  the 
theoretical  values  in  a  Monte  Carlo  experiment  using  the  screen  generator.  Monte  Carlo 
experiments  performed  using  the  screen  generator  showed  the  amplitude  perturbations 
can  significantly  reduce  the  accuracy  of  full-aperture  tilt  estimation  using  image  centroid 
motion.  However,  since  the  amplitude  perturbations  affect  the  image  centroid,  the  tilt 
estimate  using  the  combined  amplitude  and  phase  screen  allows  for  a  higher  Strehl  ratio 
than  using  only  the  phase  to  estimate  the  correction. 
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ATMOSPHERIC  TURBULENCE  SCINTILLATION  EFFECTS  ON 


WAVEFRONT  TILT  ESTIMATION 

I.  Introduction 

When  the  aperture  size  of  an  optical  system  limits  the  resolution  of  the  system  it 
is  referred  to  as  diffraction  limited.  Atmospheric  turbulence  replaces  the  aperture  size 
as  the  limiting  factor  in  optical  performance  for  large  aperture  systems,  when  imaging 
or  propagating  through  the  earth’s  atmosphere.  The  heating  and  cooling  of  the  earth’s 
surface  along  with  a  mixing  eflfect  of  wind  causes  the  atmosphere  to  be  tmrbulent.  As  a 
wavefront  from  an  object  travels  through  this  turbulent  atmosphere,  the  turbulence  will 
cause  different  sections  of  the  wavefront  to  experience  different  optical  path  lengths.  The 
optical  path  diflferences  will  bend  the  wavefront,  resulting  in  optical  field  perturbations. 
These  random  fluctuations  in  the  optical  field  can  cause  extensive  blurring  in  imaging 
through  the  atmosphere  and  is  the  limiting  factor  in  large  aperture  ground-based  imaging 
systems.  These  fluctuations  can  also  cause  significant  spreading  of  the  spot-size  of  a  laser 
propagating  long  distances  through  the  atmosphere,  significantly  reducing  the  power  on 
target.  In  the  past  few  decades,  advances  in  adaptive  optical  systems  have  made  it  feasible 
to  approach  resolutions  near  the  diffraction  limit  for  viewing  objects  beyond  the  atmosphere 
[6,8,12]. 

Many  papers  have  developed  simulation  techniques  that  incorporate  the  phase  per¬ 
turbations  caused  by  atmospheric  turbulence  [2,4,11,12,15],  but  optical  field  perturbations 
can  take  the  form  of  both  phase  and  amplitude  perturbations.  Even  if  perfect  phase  cor¬ 
rection  is  applied  by  the  adaptive  optical  system,  it  cannot  correct  for  the  amplitude 
perturbations.  Therefore,  to  accurately  evaluate  the  atmospheric  turbulence  effects,  the 
amplitude  pertmbations  should  be  considered.  This  thesis  develops  a  random  screen  gen¬ 
erator  that  simulates  both  the  phase  and  amplitude  optical  field  perturbations  due  to  weak 
turbulence. 
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The  amplitude  perturbations  eflFects  on  imaging  have  been  studied  to  some  degree, 
but  there  hasn’t  been  much  study  in  how  the  amplitude  perturbations  affect  the  phase 
measurements  required  in  adaptive  optical  systems.  This  thesis  will  study  the  effects  of 
amplitude  perturbations  on  full  aperture  tilt  estimation.  The  tilt  estimation  in  an  adaptive 
optics  system  is  used  to  remove  the  overall  phase  tilt  in  the  perturbed  wavefront. 

1.1  Atmospheric  Turbulence  Modeling 

Modeling  the  effects  of  turbulence  has  a  received  a  great  deal  of  attention.  Most  of 
the  work  on  optical  propagation  through  turbulence  is  based  upon  the  seminal  efforts  of 
Pried,  Kolmogorov,  and  Tatarski  [3,9,13].  Kolmogorov  suggested  that  over  an  extended 
propagation  path  the  index  of  refraction  differences  could  be  modeled  as  homogeneous  and 
isotropic  random  processes.  Homogeneity  implies  the  statistics  of  the  turbulent  flow  are 
independent  of  the  position  within  the  turbulent  flow.  The  isotropic  assumption  is  much 
stronger  since  it  implies  homogeneity  and  also  requires  that  the  second  and  higher  order 
statistics  of  the  turbulence  are  only  dependent  upon  the  radial  distance  between  any  two 
points  [6,12]. 

Many  papers  have  used  the  geometric  optics  calculation  to  model  atmospheric  tur¬ 
bulence  effects  on  images  taken  through  the  earth’s  turbulent  atmosphere  [12, 15].  The 
geometric  optics  approach  models  the  turbulence  effects  by  integrating  over  the  optical 
path  differences  which  simply  add  phase  to  the  propagating  wave.  The  approach  is  often 
appropriate  when  there  is  little  far-field  tmrbulence  (i.e.,  ground  to  space  propagation). 
If,  however,  the  turbulence  is  extended  through  the  entire  medium,  the  fleld  encounters 
both  phase  and  amplitude  variations  due  to  diffraction  effects  (i.e.,  propagation  within  the 
turbulence).  In  this  thesis  we  will  incorporate  the  phase  and  amplitude  effects  caused  by 
atmospheric  turbulence. 

The  phase  perturbations  have  the  effect  of  creating  a  wavefront  randomly  crinkled  in 
appearance  and  propagating  in  a  random  direction  (or  tilt).  The  amplitude  effects,  often 
referred  to  as  scintillation,  create  a  speckled  appearance  in  the  intensity  of  the  fleld.  The 
amplitude  and  phase  perturbations  have  the  effect  of  degrading  the  image  and  widening  the 
point  spread  function  (psf).  The  random  overall  tilt  of  the  wavefront  phase  perturbations 
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has  the  effect  of  moving  the  image  around  in  space,  but  does  not  reduce  its  sharpness.  If  a 
very-short-exposure  image  is  recorded  (i.e.,  the  tilt  is  removed),  the  psf  can  be  significantly 
reduced  and  resolution  enhanced.  A  large  portion  of  the  image  degradation  is  due  to 
the  random  tilt  portion  of  the  phase  perturbations  of  the  wavefront  at  the  pupil  [6,12]. 
Since  adaptive  optics  systems  can  only  correct  the  phase  perturbations,  scintillation  effects 
represent  the  theoretical  limit  one  can  achieve  in  adaptive  optics  systems. 

First,  the  phase  and  amplitude  effects  of  the  random  fluctuations  on  the  optical  field 
propagation  must  be  predicted.  The  analysis  presented  here  considers  a  diffraction  calcu¬ 
lation  based  upon  the  work  of  Lee  and  Harp  to  determine  both  the  phase  and  amplitude 
perturbation  effects  caused  by  atmospheric  turbulence  on  the  wavefront  in  the  aperture. 
The  optical  field  statistics  of  the  phase  and  amplitude  perturbations  used  in  this  thesis 
are  taken  from  the  correlation  functions  developed  by  Lee  and  Harp  [10].  In  addition,  the 
cross  correlation  of  the  phase  and  amplitude  perturbations  is  derived  for  the  first  time 
in  this  thesis.  The  correlations  are  determined  by  modeling  the  turbulence  as  layers  of 
randomly  varying  refractivity  perpendicular  to  the  propagation  path.  As  the  field  prop¬ 
agates  through  the  medium,  diffraction  occurs  at  each  of  the  layers.  Due  to  the  weak 
turbulence  assumption,  the  effects  of  the  diffracted  wave  being  diffracted  again  (second 
order  effects)  by  subsequent  layers  is  ignored.  The  weak  index  of  refraction  fluctuation 
model  does  not  predict  real  effects  of  strong  turbulence,  such  as  the  saturation  of  the  log 
amplitude  fluctuations.  To  represent  the  pertmrbations  in  the  frequency  domain  the  power 
spectral  densities  (psd)s  are  determined  by  taking  the  Fomier  transform  of  the  correlation 
functions. 

Analysis  of  the  correlation  functions  often  results  in  difficult  analytical  solutions 
with  no  closed  form.  Therefore,  simulations  of  the  atmospheric  turbulence  are  often  used 
to  determine  how  optical  systems  will  perform  [2,4,11,12,15].  These  previous  works  only 
accounted  for  the  phase  effects  of  the  turbulence.  To  investigate  the  effect  of  propagation 
through  extended  turbulence,  amplitude  effects  must  be  incorporated  as  well. 

In  order  to  study  the  effects  of  the  amplitude  and  phase  perturbations  screens  will 
be  generated  to  simulate  the  complete  random  optical  field  caused  by  the  turbulence.  A 
new  method  of  atmospheric  screen  generation  is  developed  for  use  in  performance  calcula- 
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tions  of  adaptive  optics  systems  that  is  valid  over  a  wide  range  of  atmospheric  turbulence 
parameters  and  incorporates  both  phase  and  amplitude  effects.  The  generated  random 
optical  field  over  the  aperture,  E{x),  can  be  represented  as 

E(x)  =  (l  +  i(x))eJ^(^),  (1.1) 

where  A{x)  represents  the  amplitude  perturbations  and  ^{x)  represents  the  phase  pertur¬ 
bations.  The  screen  generator  accounts  for  diffraction  effects  caused  by  turbulence  and 
incorporates  the  phase,  amplitude,  and  cross  statistics  of  the  weak  tmbulence  model.  The 
screen  generation  approach  follows  Welsh’s  paper  on  Fourier  series  (FS)  based  phase  screen 
generation  [15].  Welsh’s  relationship  for  the  phase  screen  generation  is  extended  by  using 
the  phase,  amplitude,  and  cross  psds.  A  modal  representation  is  used  to  perform  a  FS 
expansion  of  the  wavefront  phase  and  amplitude  over  a  square  area.  The  screen  generator 
uses  the  psds  of  the  phase  and  amplitude  perturbations,  determined  from  the  correlation 
functions,  in  the  FS  expansion.  The  screen  also  incorporates  the  cross  psd  derived  in  this 
thesis  to  ensiure  the  cross  statistics  are  satisfied. 

1.2  Adaptive  Optic  Systems 

Adaptive  optic  systems  operate  by  measuring  the  wavefront  phase  of  a  reference 
guide  star,  laser  guide  star,  or  the  glint  oflF  the  object  to  be  “imaged”  and  subtracting  the 
phase  from  the  object  wavefront,  resulting  in  an  improvement  in  image  quality  [12,14].  On 
many  systems,  the  wavefront  measurements  are  made  with  a  wavefront  sensor  (WFS)  such 
as  the  Shack-Hartman  sensor  and  a  separate  full-aperture  tilt  sensor.  The  measurements 
of  the  WFS  and  the  tilt  sensor  are  converted  to  electrical  signals  which  drive  actuators  on 
the  deformable  mirror  and  the  tip-tilt  mirror,  resulting  in  the  desired  phase  correction.  A 
simplified  adaptive  optics  system  is  illustrated  in  Fig.  1.1.  When  the  phase  is  perfectly 
measured  and  compensated  by  a  system  for  a  beacon  along  the  same  path  as  the  object, 
theory  predicts  perfect  phase  correction  [12].  The  result  of  perfect  correction  in  a  phase 
only  perturbed  wavefront  is  referred  to  as  diffraction  limited  performance. 

However,  adaptive  optic  systems  can  never  reach  perfect  correction.  Errors  due  to 
sampling  and  wavefront  reconstruction  will  always  exist.  In  addition,  when  the  object  and 
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Incident  Wave  Front 


Figure  1.1  Simplified  optical  configuration  of  an  adaptive  optics  imaging  system 

reference  beacon  are  not  collocated  the  light  travels  through  different  turbulence.  Since  the 
system  is  compensating  for  the  beacon’s  tmbulent  path  and  not  the  turbulence  the  object 
sees,  errors  in  the  compensation  result.  This  is  called  anisoplanatism.  Furthermore,  even 
if  perfect  phase  correction  is  applied,  the  adaptive  optics  system  cannot  correct  for  the 
amplitude  perturbations.  Therefore,  the  theoretical  limit  of  adaptive  optic  compensation 
is  determined  by  these  scintillation  effects.  To  accurately  evaluate  total  effects  of  the 
atmospheric  turbulence  the  scintillation  efiects  should  be  considered. 

Scintillation  effects  in  imaging  have  been  studied  to  some  degree,  but  there  hasn’t 
been  much  study  in  how  the  amplitude  perturbations  affect  the  phase  measurements  re¬ 
quired  in  adaptive  optic  systems.  This  thesis  studies  the  effects  of  scintillation  on  full 
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aperture  tilt  estimation  in  weak  turbulence.  The  tilt  estimation  is  used  by  the  tip-tilt  mir¬ 
ror  in  an  adaptive  optics  system  to  remove  the  overall  tilt  in  the  wavefront.  Monte-Carlo 
experiments  using  the  random  screen  generator  will  be  performed  to  investigate  the  effects 
of  the  amplitude  perturbations  on  the  full-aperture  tilt  estimation.  Since  the  random  tilt 
contains  87%  of  the  power  in  the  wavefront  phase  fluctuations,  accurate  wavefront  tilt  es¬ 
timation  is  essential  to  an  adaptive  optics  system.  These  experiments  will  use  the  random 
amplitude  and  phase  perturbation  screens  produced  by  the  screen  generator  developed  in 
this  thesis  to  determine  the  effect  of  the  amplitude  pertmbations  on  the  tilt  estimation. 

1.3  Summary  of  key  results 

The  theoretical  statistics  were  developed  for  the  amplitude  and  phase  perturbations 
of  the  optical  field  caused  by  weak  turbulence.  The  statistics  were  applied  in  the  develop¬ 
ment  of  a  random  optical  field  screen  generator.  The  mean  square  value  and  the  structure 
functions  of  the  phase,  amplitude,  and  cross  terms  are  calculated  in  a  Monte  Carlo  exper¬ 
iment  using  the  screen  generator  to  within  1%  of  the  theoretical  values  for  a  number  of 
screen  dimensions.  Monte  Carlo  experiments  performed  using  the  screen  generator  deter¬ 
mined  the  amplitude  perturbations  can  significantly  reduce  the  accuracy  of  full-aperture 
tilt  estimation  using  image  centroid  motion.  However,  using  the  combined  amplitude  and 
phase  screen  image  centroid  to  estimate  tilt  correction  allowed  for  a  higher  Strehl  ratio 
than  using  only  the  phase  to  estimate  the  correction. 

1.4  Overview 

Chapter  II  develops  the  amplitude,  phase,  and  cross  correlation  functions  of  the  at¬ 
mospheric  turbulence.  Chapter  III  determines  the  power  spectral  density  (psd)  functions 
of  the  atmospheric  turbulence  from  the  correlation  functions.  Chapter  IV  uses  the  psds 
developed  in  Chapter  III  to  develop  an  atmospheric  turbulence  phase  and  amplitude  screen 
generator.  Chapter  V  studies  the  effect  of  the  amplitude  perturbations  on  the  wavefront 
tilt  estimation  and  the  resulting  Strehl  ratios.  Results  are  presented  that  indicate  how  am¬ 
plitude  perturbations  affect  the  tilt  estimation  for  various  conditions  of  weak  atmospheric 
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turbulence.  Chapter  VI  gives  some  conclusions  about  the  results  and  recommendations  for 
further  study. 
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11.  Correlation  Functions  of  Weak  Turbulence 

The  purpose  of  this  section  is  to  develop  the  phase,  amplitude  and  cross  correlation  func¬ 
tions  for  the  optical  field  perturbations  caused  by  atmospheric  turbulence.  A  diffraction 
based  calculation  is  used  to  develop  the  correlation  functions  associated  with  weak  scat¬ 
tering  through  a  random  medium.  The  approach  follows  Lee  and  Harp’s  development  and 
results  in  the  same  functions  as  in  their  paper  [10].  In  addition  to  the  phase  and  amplitude 
correlations,  the  cross  correlation  is  derived  for  the  first  time  and  follows  an  approach 
similar  to  Troxel’s  paper  [14]. 

The  turbulence  induced  random  index  of  refraction  is  modeled  with  thin  slabs  of 
random  refractivity  perpendicular  to  the  propagation  path.  Each  slab  is  represented  with 
a  Fourier  decomposition  of  random  phase  and  varying  spatial  frequency  (see  Fig.  2.1).  A 
wave  is  propagated  from  slab  to  slab  where  each  slab  can  be  treated  as  a  phase-diffraction 
grating.  At  each  slab,  the  incident  wave  is  diffracted  resulting  in  two  diffracted  orders 
and  the  original  direction  of  propagation.  Due  to  a  weak  scattering  approximation  only  a 
single  diffraction  at  each  slab  is  allowed  (i.e.,  the  diffracted  orders  aren’t  allowed  to  diffract 
in  subsequent  layers).  Therefore  only  the  zeroth  order  has  enough  power  to  have  a  non- 
negligible  effect  at  the  next  layer.  The  results  are  integrated  over  the  spatial  frequencies 
and  every  angle  of  the  Fourier  components  to  determine  the  amplitude  and  phase  pertur¬ 
bations  [10].  Finally,  the  random  field  perturbations  in  the  amplitude  and  phase  are  used 
to  calculate  the  amplitude,  phase,  and  cross  correlations. 

The  final  solutions  for  the  amplitude,  phase  and  cross  correlations  are  very  similar 
in  form  and  each  include  three  multiplicative  terms:  a  term  relating  the  perturbations  at 
one  point  with  those  at  another  point,  the  power  spectrum  of  the  refractive  index,  and  a 
weighting  term  which  accounts  for  the  relative  perturbations  at  a  given  position  along  the 
transmission  path  [10]. 

2.1  Derivation  of  Correlation  Functions 

Using  the  layered  approach,  the  effect  on  the  optical  field  can  be  represented  as  a 
differential  amplitude  perturbation  and  a  differential  residual  phase  at  each  slab  or  layer 
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Figure  2.1  A  cross  section  of  the  Fourier  representation  of  the  random  refractivity  slab, 
where  A{n,  k)  is  the  amplitude,  (poiVy  is  the  random  phase  and  ^  is  the 
spatial  frequency. 

(see  Fig.  2.1).  The  total  amplitude  pertmrbation  represents  the  deviation  from  a  unit 
amplitude.  The  total  phase  perturbation  represents  the  total  phase  deviation  induced  by 
the  refractive  index  perturbation.  Lee  and  Harp  derives  the  following  differential  phase 
and  amplitude  perturbations  for  propagation  through  weak  turbulence  modeled  as  the 
propagation  through  one  of  the  differential  slabs  [10] 

dPa  =  —ka{K)dr]  cos{k{x  +  b))  —  fc)]  (2.1) 

and 

dPp  =  —ka{K)d7]  cos(K(a;  +  b))  cos[i]{Vk^  —  —  ^)].  (2.2) 

Lee  and  Harp  define  k  as  the  wavenumber,  27r/A;  a(K)  as  the  peak  amplitude  of  the 
refractive  index  perturbations;  b  as  the  random  phase  of  the  refractive  index;  drj  as  the 
differential  propagation  path  length;  rj  as  the  path  length  from  the  slab  to  the  observation 
plane;  x  as  the  direction  along  which  the  periodic  Fourier  component  is  taken;  and  k  as 
the  spatial  frequency  of  the  particular  Fourier  component. 
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Using  the  binomial  expansion  and  the  assumption  that  k  -C  the  following  factor 
can  be  simplified  to 


y/k"^  —  —  k  =  k{J  1  —  ^  —  1)  «  fc 


K 


/  k^l  ^ 


2k 


a. 


(2.3) 


Now,  Troxel  represents  the  perturbations  in  three  dimensions  at  a  particular  altitude,  t/, 

by 


Paiv,  K)  x)  =  A{t],  k)  cos(/c  •  X  +  (I>o{t),  k))  sm{T]a) 


(2.4) 


and 


Pp{ri,  K,  x)  =  —A(ri,  k)  cos(k  ■  x  +  <j>o{Vi «))  003(770;),  (2.5) 

where  A(77,  /c)  and  <i>oiVj  k)  are  related  to  the  amplitude  and  phase  of  the  refractivity  field 
Fourier  component  and  x  is  a  position  vector  in  the  plane  of  the  slab  [14].  The  total  residual 
phase  and  amplitude  perturbations  are  found  by  integrating  over  all  slabs  distributed  along 
a  path  length,  L,  and  all  frequencies,  it  and  can  be  represented  as 

rL 


Pai^)  ~  j  j  d77dKA(77,  k)  cos(-k  •  x  —  ^o(’?j «))  sin(77a) 


(2.6) 


Pp{x)  =  —  J  J  dT]dKA{T],K)cos{—K’X  —  <j>o{r],K))cos{r]a),  (2.7) 

where  o  =  K?/{2k).  Representing  the  perturbations  in  complex  phasor  notation,  the 
amplitude  perturbations  and  residual  phase  perturbations  of  the  field  can  be  written  as 


Pa{^)  —  —P'{  I  I  dr]dKA{T],K)e 


p-j  («•?) 


smi 


(77a) 


and 


Pp{x)  = j  dr]dKA{r],K)e  cos(77a)| , 


(2.8) 


(2.9) 


where  K  is  an  operator  that  takes  the  real  portion  of  a  complex  quantity.  The  TZ  operator 
is  now  dropped  and  the  work  continues  with  the  complex  phasor  representations  Pa{x)  and 
Ppjx)  given  by 


Pa{x)  —  J  J  d7)ditA{Tj,  K)e  sin(77Q:) 


(2.10) 
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and 


Pp{x)  =  j  j  dr}dKA{T],K)e~^^^'^  cos{r]a),  (2-11) 

where  A{ri,  it)  =  A{ri,  and  a  =  K^/{2k)  [14]. 

Prom  Eqs.  (2.10)  and  (2.11)  above  the  development  of  the  amplitude  and  phase 
correlations  directly  follows  Troxel’s  approach  [14].  The  amplitude  and  phase  correlations 
are  defined  as  the  ensemble  average  of  the  perturbations  at  two  different  locations  separated 
by  d  and  can  represented  as 


raa(^  =  {Pa(x)Pa(x-^) 

and 


(2.12) 


Tpp{d)  =  (Pp{S)Pp{x  -  dl)  ,  (2.13) 

where  the  angle  brackets  denote  the  expected  value  of  the  argument. 

Lee  and  Harp  determines  the  correlations  of  the  amplitude  and  phase  perturbations 
in  weak  turbulence  as  the  following 

rQa(d)  =  4fc^7r^  /  /  d/cd7j/cJo(/cd)$(K)sin^(Q;??),  (2.14) 

Jo  Jo 

and 


rpp(d)  =  4fc^7r^  /  f  dKdTjKJo{Kd)^{K)cos^{ar)), 
Jo  Jo 


(2.15) 


fL  roo 

to  Jo 

where  a  =  fjg  and  Jq  is  the  zeroth  order  Bessel  function  of  the  first  kind  [10].  Equations 
(2.14)  and  (2.15)  are  the  key  equations  that  describe  the  second  order  statistics  of  the 
amplitude  and  phase  effects  due  to  weak  turbulence.  These  equations  form  the  basis  for 
all  other  statistical  representations  of  the  field  effects  due  to  the  atmospheric  tmrbulence 
in  this  thesis. 


Another  key  second  order  statistic  is  required  to  define  the  atmospheric  turbulence 
effects  on  the  optical  field.  In  order  to  investigate  how  the  phase  and  amplitude  perturba¬ 
tions  are  related  one  needs  to  know  how  they  are  correlated  (i.e.,  their  cross  correlation). 
The  cross  correlation  is  defined  as  the  ensemble  average  of  the  amplitude  and  phase  per¬ 
turbations  taken  at  two  different  locations  separated  by  d.  The  method  used  here  is  very 
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similar  to  the  approach  used  by  Troxel  to  derive  the  angle  dependent  phase  and  amplitude 
correlation  functions  [14].  Using  the  perturbation  relationships  in  Eqs.  (2.10)  and  (2.11), 
the  cross  correlation  can  be  determined  by 


TaM  =  {Pa(x)P;ix-d)), 

=  (^J  j  Jq  d??ld7?2dKld^^A(7?l,/cl)e“^’(''l■‘^sin(77la:l) 

xA*{t]2,  cos(r)20t2)^ . 


(2.16) 


Since  the  only  random  quantities  axe  in  the  A  terms,  the  other  deterministic  terms 
can  be  passed  through  the  expectation. 


^ap{^  =  j  j  Jq  <i’?l<l’?2dKldK2^i(77l,Kl)A*(7?2,/C2)) 

cos{r}2a2)- 


(2.17) 


Troxel  uses  the  following  relationship  derived  from  a  Tatarski  relation  to  perform  the 
expectation: 


^i(7?i,Ki)i*(772,K2))  =  87rA:^5(Ki  -  K2)S{vi  -  ?72)^(«1,’7i), 


(2.18) 


where  $(«,  77)  is  the  three  dimensional  refractive  index  power  spectrum  with  +  Ky, 

and  6  is  the  Dirac  delta  function  [13, 14].  The  ^(771  —  772)  part  of  Eq.  (2.18)  is  a  derived 
relation  and  should  not  be  considered  as  a  statement  of  independence  of  turbulent  layers. 
Substituting  Eq.  (2.18)  into  Eq.  (2.17)  gives 


Tap{^  =  Sk'^TT  j  J  dTjdK^{K,ri)sm{ar})cos{a'i])e 

and  using  a  trigonometric  identity  the  relation  becomes 

Tapi^  =  4:k^‘iv  J  J  dr]dK^{K,r])sm{2aT])e~^^'^. 


(2.19) 


(2.20) 
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The  following  relation  from  Goodman  illustrates  that  the  real  part  of  the  complex  corre¬ 
lation  is  twice  the  real  valued  cross  correlation: 


r,..(T)  =  2r(:/)(r)  -I- (2.21) 

where  is  the  real  correlation  [6].  Therefore,  the  real  correlation  is  half  the  real 

portion  of  the  complex  one  (i.e.  Tap{^  =  5^{rap(d)}).  The  real  correlation  can  now  be 
represented  as 


Tap(^  =  2k^'K  j  J  d7?dK$(/c,77)sin(2a7/)cos(K  •  d). 


(2.22) 


The  vector  representation  can  be  changed  to  polar  coordinates  where  dw  =  «:d/cd0. 
The  integration  over  6.6  goes  from  0  to  tt.  The  integration  over  69  can  then  be  performed 
using  the  relation 


d0cos(acos(0))  =  7rJo(a),  (2.23) 

where  Jq  is  the  zeroth  order  Bessel  function  of  the  first  kind.  Finally,  the  cross  correlation 
can  be  represented  as 


rap(<^  =  2A:^7r  f  f  f  dfld/cdTj/ccos  (Kdcos(0))  $(k,  77)  sin(2Q;77) 
Jo  Jo  Jo 

rL  roo 

=  2fc^7r^  /  /  6K6T]KJo{Kd)^{K,T])sm{2ar]). 

Jo  Jo 


(2.24) 


Equation  (2.24)  is  the  cross-correlation  of  the  amplitude  and  phase  perturbations  of  the 
optical  field  due  to  weak  turbulence.  This  equation  will  be  used  throughout  this  thesis  to 
represent  statistics  of  the  cross  terms  of  the  phase  and  amplitude. 

2.2  Correlations  for  a  constant  strength  of  turbulence 

In  order  to  proceed  with  the  correlations  developed  in  this  section  they  need  to 
be  related  to  a  particular  turbulent  medium.  The  3-D  turbulence  power  spectrum  of 
the  refractive  index,  $(/c,  77),  is  a  function  of  both  the  turbulence  along  the  path  and 
the  frequency  distribution.  The  following  relation  essentially  breaks  up  the  functional 
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dependence  of  the  3-D  power  spectrum  into  the  turbulence  strength  profile  and  a 

normalized  power  spectral  density  can  be  represented  as  the  following  [6,12] 

$(«,77)  =  0.033C'2(,7)$o(«),  (2.25) 

where  represents  the  strength  of  turbulence  along  the  path.  The  von  Karman 

turbulence  power  spectrum  definition  for  ^o(«)  uses  an  inner  scale,  Lm,  and  outer  scale, 
Lo,  to  represent  the  smallest  and  largest  turbulence  scale  sizes,  respectively,  in  a  particular 
medium.  Lg  is  usually  on  the  order  of  ten  to  twenty  meters,  where  as  Lm  is  usually  only 
a  few  millimeters  [12].  The  von  Karman  spectrum,  $o(k),  is  defined  as  [12] 


2  t2 

e 


(k2  +  ^)11/6' 


(2.26) 


Using  Eq.  (2.25),  the  amplitude,  phase  and  cross  correlations  given  in  Eqs.  (2.14), 
(2.15),  and  (2.24)  can  be  written  as 


nL 

dj;dKK Jo(/cd)$(K,  ri)  sin^(Q;77) 

=  4/i:^7r^ (0.033)  J  dKKjQ{Kd)^o{K)  j  dr)C^{r})six^{ar]), 


(2.27) 


and 


njj 

dTjd/cK  Jo(Kd)$(«;,  if)  cos^  {at]) 

_ 

=  4fc^7r^(0.033)  J  dK«Jo(Kd)$o(K)  J  dT]C^{T))  cos^ {ar]) , 


(2.28) 


and 


rap(<^  =  J  J  dT]dKKJo{Kd)^{K,r])sm{2aT)) 

roo  rL 

=  2A:^7r^ (0.033)  /  dKKJo(Kd)$o(«)  /  d77C^(77)  sin(2Q;?7). 
Jo  Jo 


(2.29) 


To  compare  the  correlations  above  let’s  assume  the  strength  of  turbulence  is  constant 
over  the  entire  length  of  propagation.  For  horizontal  propagation  through  the  atmosphere 
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constant  turbulence  is  a  good  approximation,  but  for  propagation  up  through  the  atmo¬ 
sphere  it  is  not  such  a  good  approximation,  since  strength  of  turbulence  has  a  strong 
dependence  on  altitude.  The  total  turbulence  over  the  entire  path  is  defined  as 

(2-30) 

where  tq  is  the  Pried  parameter  of  the  turbulence  and  A  is  the  average  wavelength  of  the 
light  [3].  The  Pried  parameter  can  be  interpreted  as  the  limiting  aperture  size  beyond 
which  no  increase  in  resolution  occurs. 

Using  Eq.  (2.30),  now  the  correlations  can  be  integrated  over  the  path  length  and 
calculated  as  the  following 


\2  roo  rL 

^aai^  —  (0.033)— ( -^ - )  /  d/CKJo(/cd)$o(«)  /  dr}sm^{ar]) 

L  \  tq  J  Jo  Jo 

/niRfiNS/s  roo  rr  1 

=  4A:^7r^(0.033)-— (-^ - )  /  dK/cJo(/cd)$o(K)  —  —  —  sin(2Q'L) 

L  \  fQ  J  Jo  [2  4q! 

(0 185  \  r  1  1 

- )  /  dKKJo(Kd)$o(^)  1  — )  (2.31) 

ro  J  Jo  L  2q'L  J 


and 


\2  /Q  18C>\5/3  i>oo  rL 

rpp(d)  =  4A:^7r^  (0.033)—  - j  J  d«;/cJo(/cd)$o('«)  J  drj  cos^  {arj) 

\2  foo  T  L  ^ 

=  4A:^7r^(0.033)— ( - )  /  dKKJo('^®^)^o('<^)  x  +  sin(2aT) 

J^  \  ro  /  Jo  L2  4q; 


/0  185\®/^  roo  r  1 

=  Stt"* (0.033)  ( —^ - )  /  dK«Jo(red)$o('^)  1  +  7; — ;r  sin(2aL) 

\  ro  J  Jo  L  2oiL 


,  (2.32) 


and 


no  /0  185\®/^  r°° 

T'apW  —  2A:^7r^(0.033)-jr  I - )  /  dK/cJo(/cd)$o(K)  /  d?; sin(2Q:7?) 

\  ro  J  Jo  Jo 

\2  /ni85\V3  roo  1 

=  2fc^7r^(0.033)—  ( -^ - )  /  dKKJo(Kd)$o(K)~  sin^(a:Z’) 

\  ro  /  Jo  ot 

/O  185  \  r°°  1 

=  87r'^(0.033)  I - )  /  dKKJo{Kd)^o{i^)—rsin^{aL),  (2.33) 

\  ro  /  Jo  (xL 
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where  o  =  If  =  ^. 

In  order  to  present  these  correlations  more  clearly,  it  is  also  convenient  to  perform 
a  change  of  variables  when  calculating  the  correlations.  Let  k'  =  kLq.  Performing  the 
change  of  variables  results  in 


27rL2  .  /aLk'2\ 


and 


(T  \  5/3  roo  1 

(K'd\  r  27rLo  .  f  XLk'^W 


and 


r„p(d)  =  87r'^(0.033)(0.185)®/3 


47r2)ii/6 


(2.34) 


(2.35) 


(2.36) 


where  we  have  approximated  the  exponential  in  Eq.  (2.26)  as  being  approximately  one. 
Numerical  integration  was  performed  over  k'  for  the  correlations  above  for  a  strength  of 
turbulence,  Lo/rg  =  100,  and  ^  =  20000.  The  correlation  plots  are  shown  in  Figs.  2.2, 
2.3,  and  2.4.  Notice  how  much  wider  rpp(<l)  is  than  roo(<^  and  Tap{^-  This  trend  indicates 
that  the  phase  perturbations  are  spatially  correlated  over  a  larger  area  than  the  amplitude 
perturbations. 
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Figure  2.2  The  normalized  amplitude  correlation  is  plotted  versus  for  a  strength  of 
turbulence,  Lofro  =  100,  and  ^  =  20000. 


Figure  2.3  The  normalized  phase  correlation  is  plotted  versus  for  a  strength  of  tur¬ 
bulence,  Lo/rg  =  100,  and  ^  =  20000. 
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2.3  Structure  Functions 


The  structure  functions  for  these  perturbations  can  be  represented  as  [6, 12] 


Daai^  =  2Taa(0)-2Taai^, 

=  16,‘(0.033)(0.185)V3  ^  + 


and 


and 


-Me 


XLk^ 


Sin 


(XLk'^\ 

I  2^Ll  I 


Dpp(i)  =  2rpp(o)  -  2rpp(^, 

-  '“‘"“■iww"  &rT 


,  27rL2  .  /ALk'2\ 

^  ^  AL«'2  27rT2  j  J  ’ 


(2.37) 


(2.38) 


Dap{^  —  2r(jp(o)  —  2rop(<i), 

/  T  \  V®  roo  1 

=  16.‘(0.033)(0.185)V3  (^)  ^ 


X 


I- ‘©I 


(SS)-'— 


47rL2 


(2.39) 


The  wave  structure  function  is  the  sum  of  the  phase  and  amplitude  structure  functions 

[6,12] 


D[S)  —  Dpp(Jt)  +  Daa{^ 

=  32.3(0.033)(0.185)3/3 
L  ,  /K'dM 


(2.40) 
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2.4  Mean  Square  Value  Functions 

The  mean  square  value  of  the  amplitude  and  phase  can  be  calculated  easily  from  the 
correlations.  Letting  d  go  to  zero  in  Eqs.  (2.34),  (2.35),  and  (2.36)  we  get  the  ensemble 
average  of  the  squared  phase  and  amplitude  perturbations  or  the  mean  square  values. 
These  can  be  represented  as  follows 


TaaiO)  =  87r4(0.033)(0.185)5/3 


47r2)ii/6 


and 


and 


,  27rL2  .  (\Lk!^\ 


r'pp(o)  — 


,  /  r  \  5/3  roo  1 

8.‘(0.033)(0.185)'/3  ^ 

,  2t: Ll  .  /ALk'2\] 


/r  \  too  1 

r.,(0)  =  8^‘{0.033)(0.185)=/= 

/47rL2\  .  2  /AL/c'2\ 

^  (ALK'2j®i“  )' 


(2.41) 


(2.42) 


(2.43) 


Normally,  as  the  propagation  path  or  wavelength  increases  Lo/vo  increases  and  the 
turbulence  effects  are  stronger.  For  the  following  plots  the  strength  of  turbulence,  Lo/^oj 
is  held  constant  to  show  the  relative  mean  square  values  of  the  phase  and  amplitude 
perturbations.  Figure  2.5  shows  the  mean  square  value  of  the  amplitude  perturbations  as 
a  function  of  ^  for  different  strengths  of  tmbulence,  Lo/to.  For  constant  Lolvo,  notice 
how  the  amplitude  perturbations  decrease  substantially  as  gets  large.  On  the  other 

r  2 

hand,  the  phase  contribution,  plotted  in  Fig.  2.6,  gets  very  large  for  large  Note  also, 

r  2 

the  sum  of  the  phase  and  amplitude  mean  square  values  is  invariant  to  for  constant 
Lo/tq.  For  instance,  for  Lo/rg  =  100  the  sum  of  the  phase  and  amplitude  mean  square 
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values  is  186.65.  The  mean  cross  value,  rap(O),  of  the  phase  and  amplitude  perturbations 

L/^  * 

is  plotted  in  Fig.  2.7.  Notice  how  rap(O)  decreases  as  ^  increases  similar  to  the  amplitude 
perturbations  in  roo(O). 


Figure  2.5  The  mean  square  value  of  the  amplitude  perturbations  in  the  pupil  are  plotted 
for  different  strengths  of  turbulence,  Lo/r,,. 
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Figure  2.6  The  mean  square  value  of  the  phase  perturbations  in  the  pupil  are  plotted 
for  different  strengths  of  turbulence,  Lo/ro- 


XL 

Figure  2.7  The  mean  cross  value  of  the  phase  and  amplitude  perturbations  in  the  pupil 
are  plotted  for  different  strengths  of  tiurbulence,  Tq/^o- 


III.  Phase,  Amplitude,  and  Cross  Power  Spectral  Densities  of  Weak  Turbulence 


Often  it  is  convenient  to  represent  the  perturbations  in  the  frequency  domain  to  get  a  better 
idea  of  the  strength  of  tmbulence  and  the  perturbations  at  the  pupil.  The  frequency  domain 
representations,  the  power  spectral  densities  (psds),  are  developed  from  the  correlations 
introduced  in  the  previous  section.  The  psds  are  calculated  by  taking  the  Fourier  transform 
of  the  correlation  functions  [6].  The  psds  developed  in  this  chapter  will  be  used  to  create 
random  phase  and  amplitude  field  screens  with  the  correct  first  and  second  order  statistics 
in  the  next  chapter. 


To  begin,  below  are  the  phase,  amplitude  and  cross  correlation  functions  of  the  field 
perturbations  due  to  the  refractive  index  perturbations  that  were  developed  in  Eqs.(2.15), 
(2.14),  and  (2.24)  in  the  last  chapter 


and 


and 


rL  roo 

Uaa{.d)=A.k‘^'K^  1  /  d«d7;/cJo(/cd)$(/c)sin^(a»7), 

Jo  Jo 

(3.1) 

Tpp{d)  =  f  [  dKdr]KJo{Kd)^{K)cos^{ar)), 

Jo  Jo 

(3.2) 

Tap{d)  =  2k^'K^  J  J  dKd77KJo(Kd)#(K)  sin(2a77). 

(3.3) 

11 

The  definitions  of  the  Fourier  Bessel  transform  and  inverse  transform  are  the  follow- 


ing  [6| 


m  =  ^  {/(»;)} 

oo 

27rda;/  {x)xJq  {2Trpx) 

and 


/w 


{F(p)} 

roo 

/  2TrdpF{p)pjQ{2Trpx). 

Jo 


(3.4) 


(3.5) 
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The  Fourier  Bessel  transform  is  used  to  perform  a  Fourier  transform  of  a  function 
represented  in  polar  coordinates.  If  the  correlations  are  represented  in  polar  coordinates, 
the  Fourier  Bessel  transform  definition  above  can  be  utilized  to  perform  the  Fourier  trans¬ 
form.  After  making  a  change  in  variables  to  polar  coordinates,  where  k  =  27rp,  d/c  =  27rdp, 
and  a  =  (47r^  |  p  p)/(2A:),  the  correlation  function  for  the  amplitude  can  be  written  as  the 
following 


rL  foo 

^aa{d)  =  4fc^7r^  dri  2'Kdp^{2'!rp,r))sm^{a'q)2iTpJo{2'Kpd) 
Jo  Jo 

=  J  |$(27rp,?j)sin^(Q:77)| . 


(3.6) 


Utilizing  the  Wiener  Khinchin  theorem,  the  psd  can  be  found  by  taking  the  Fourier 
transform  of  the  correlation  function  [6].  Performing  the  Fourier  transform  on  the  ampli¬ 
tude  correlation  function  above  yields  the  following  equation: 


$aa(p)  =  H^ppid)} 

=  ^  J  |$(K,??)sin^(a7?)|| 

fL 

=  8k^7T^  /  dr]^{2np,r])sm^(aij). 

Jo 


(3.7) 


After  making  the  same  change  in  variables  as  before,  the  phase  correlation  can  be 
written  as  the  following 


rL  poo 

rj,j,(d)  =  Ak^TT^  dr]  2Trdp^{2Trp,T])cos^{aTj)2'KpJo{27rp  d) 

Jo  Jo 

=  Sk^TT^  J  d7?.F~^  |$(27rp,7/)cos^(Q;77)| , 

where  re  =  2ttp  and  dre  =  27rdp  and  now  a  =  (47r^  |  p  p)/(2fc). 


(3.8) 
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Similarly,  applying  the  Fourier  transform  to  the  phase  correlation  function  above 
yields  the  following  equation  for  the  phase  psd 


^pp{p)  =  •^{r'pp(d)} 

=  J  ■!$(«,  7/)  cos^(a7?)j 

rL 

=  /  dj?$(27rp,7?)cos^(Q:7/). 

Jo 


(3.9) 


Similarly,  the  cross  correlation  function  can  be  written  as  the  following 


/•L  /"oo 

FapCd)  =  2fc^7r^  /  d??  /  27rdp$(27rp,77)sin(2Q!7?)27rpJo(27rpd) 

fL 

=  4A:^7r^  /  drjJ^  ^  {$(27rp,»/)  sin(2Q!7;)}  ,  (3.10) 

Jo 


where  k  =  27rp  and  d/c  =  27rdp  and  now  a  =  (47r^  |  p  p)/(2fc).  As  before,  applying  the 
Fourier  transform  to  the  cross  correlation  function  above  yields  the  following  equation  for 
the  cross  psd: 


^ap{p)  —  •^{r'op(d)} 

=  T  i^Ak^TT^  J  dTjT~^  {^{2'Kp,r))sm{2ar])} 

rL 

=  Ak^ir^  /  d7;$(27rp,7/)sin(2a7;). 

Jo 


(3.11) 


In  all  of  the  psds  derived  above  the  3-D  psd  for  the  refractive  index  perturbations  is 


$(/c,77)  =  0.033C2(,,)$o(«), 


(3.12) 


where 


2r2 


$o(k) 


(k2  +  ^)11/6- 


(3.13) 
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Letting  k  =  2Trp  in  the  above  equation  and  approximating  the  exponential  to  one,  yields 
the  following  relationship 

(27r)-ll/3 


^o(p)  = 


(3.14) 


We  can  use  the  relationship  introduced  in  the  previous  section  for  constant  C^ir])  over  the 
propagation  path  to  represent  the  3-D  psd  for  the  refractive  index  perturbations  as 
A2 


»(„,,)  =0.033^  (2^)“'  *„(«). 

Therefore,  the  amplitude,  phase,  and  cross  power  spectral  densities  become 


^oa(/3)  = 


i-L 

8k^ir^  /  d77$(27rp,7;)sin^(Q!77) 

Jo 

fL 

8A:^7r^(0.033)$o(p)  /  ^V^niv) i^v) 

8*:V*(0.033)  f  “  i  (^)  sM2ai) 

16,r'(0.033)  (5^)  ^  *„(rt  [l  -  (2;;^)  sin(2VAL) 


and 


^pp{p)  — 


j-L 

Sk^TT^  /  d7;$(27rp,  77)008^(077) 

Jo 

rL 

8/:^7r^(0.033)$o(p)  /  d77C^ (77)  008^(077) 
^0 


=  8)fcV^  (0.033) 


A2  /0.185\5/3]  ^  ,  .\L  1  /1\  . 

(— )  J  *„W  [2  +  I  {-) 


and 


j-L 

^ap{p)  =  4A:^7r®  /  d77$(27rp,  77)  sin(2a77) 

Jo 

rL 

=  4A:^7r®(0.033)$o(/3)  /  d77C^(77) sin(2a77) 
Jo 


(3.15) 


(3.16) 


(0  185  \  r  /  1  \ 

j  $o(p)  [1  +  j  sin(27rp2AL)J  ,  (3.17) 
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=  4fc^7r®(.033)  ^  ^o(io)  ^(1  -  cos{2aL)) 

(n  1  Qc  \  5/3  r  i 

$o{p)  ^p2xi  , 


(3.18) 


where  a  =  \  p\‘^)/{2k)  and  $o(p)  =  j2+V|u/6- 

To  present  the  psds  more  clearly  a  change  of  variables  is  performed.  This  change 
of  variables  allows  the  psds  to  plotted  versus  the  dimensionless  quantity,  p'  =  pLo-  The 
$o(p0  t)e  represented  as  the  following 


/OTrl-ii/a 

but  we  find  it  convenient  to  represent  it  as  follows 

,  (2,r)-»/3i2£5/3 

*“('■)=  (/i  + 1)11/6  ■  (*■“) 

Performing  the  change  in  variables  results  in  the  following  relations  for  the  amplitude, 
phase,  and  cross  psds 


(3.20) 


^aa{p')  =  167r®(0.033)(0.185)®/3 


Lo\^/^  (27r)-iV3 

To)  (p'2  + 1)11/6 


(  Ll  \  .  /27rp'2AL 


(3.21) 


=  16’r=(0.033)(0.185)>/»(^)'‘''J|^^^ 


\  (  Ll  \  .  (2‘kp'^XL 

^  ^+(27rp'2AL)®‘“(  Ll 


(3.22) 


^ap(p')  = 


16*=(0.033)(0.185)'/’  (^)'^'  +'i)n/6 

[  -^o  •_2  ( T^p'^^LW 

i  jj' 


V  : 

7rp'2AL®'“ 


(3.23) 
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Figure  3.1  shows  the  power  spectral  densities  as  a  function  of  pLg.  Notice  how  the 
phase  is  much  stronger  than  the  amplitude  for  smaller  spatial  frequencies,  but  as  pLg  in¬ 
creases  the  phase  and  amplitude  contributions  converge  to  the  same  relationship.  Also, 
notice  how  much  quicker  the  phase  psd  and  even  the  cross  psd  drops  off  as  a  function  of  fre¬ 
quency  than  the  amplitude  psd.  The  wider  psd  of  the  amplitude  perturbations  corresponds 
to  higher  spatial  frequency  content  than  the  phase  perturbations. 


Figure  3.1  The  power  spectral  densities  as  a  function  of  pLg  for  X/Lo  =  5  x  10 
L/Lo  =  10000,  and  Lg/ro  =  100. 
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JV.  Fourier  Series  Based  Atmospheric  Phase 
and  Amplitude  Field  Screen  Generator 

Simulations  are  often  used  to  study  the  effects  of  atmospheric  turbulence  on  imaging  sys¬ 
tems  and  laser  propagation.  In  this  chapter  we  use  the  power  spectral  densities  calculated 
in  the  previous  chapter  to  develop  phase  and  amplitude  screens.  Many  screens  only  ac¬ 
count  for  phase  by  performing  a  geometric  optics  phase  calculation  [2,4, 11, 12, 15]  which 
is  appropriate  when  there  is  little  far-field  turbulence.  When,  however,  the  turbulence  is 
extended  through  the  entire  medium  the  field  encounters  phase  and  amplitude  variation 
due  to  diffraction  effects.  Using  the  proper  ensemble  statistics  of  the  amplitude,  phase, 
and  cross  power  spectral  densities  (psds)  will  allow  for  realistic  simulations  of  extended 
turbulence.  These  screens  will  simulate  wavefiront  phase  and  amplitude  perturbations  in 
the  optical  field  induced  by  propagation  through  the  turbulence. 

Many  methods  have  been  developed  to  generate  phase  screens.  The  methods  can 
be  broken  into  two  mathematical  approaches.  The  first  approach  uses  a  2-Dimensional 
rectangular  grid  of  points  to  represent  the  statistics  of  the  screen,  sample  based  screen 
representation.  Whereas,  in  the  second  approach  the  screen  is  represented  as  a  sum  of 
orthogonal  basis  functions  or  a  modal  based  representation.  The  most  common  method 
uses  an  approach  based  upon  the  Fourier  transform.  The  technique  is  direct  and  compu¬ 
tationally  efiicient  due  to  the  use  of  the  Fast  Fourier  Transform  (FFT).  The  problem  with 
this  method  is  that  the  energy  of  low  spatial  frequencies  (e.g.,  tilt)  is  under-represented. 
The  low  spatial  frequencies  contain  a  large  majority  of  the  power  which  presents  a  major 
problem  when  studying  the  effects  of  tilt.  A  model  based  approach  allows  for  much  better 
low  frequency  representation  [15]. 

The  method  used  in  this  thesis  follows  Welsh’s  paper  on  Fourier  series  based  phase 
screen  generation  [15].  A  modal  representation  is  used  to  perform  a  Fourier  Series  (FS) 
expansion  of  the  wavefront  phase  and  amplitude  over  a  square  area.  This  approach  has 
a  much  more  accurate  representation  for  low  spatial  frequencies  than  the  sample  based 
approach  [15].  Welsh’s  relationship  for  the  phase  screen  generation  is  extended  by  using 
the  phase,  amplitude,  and  cross  psds  developed  in  the  previous  chapter. 
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4.1  The  Development  of  the  Random  Screen 

To  begin,  the  complex  phase  in  the  pupil  at  the  point  x  can  be  expressed  as  an 
integral  of  the  contributions  of  each  layer  along  the  path  as  in  Eq.  (2.11): 

Pp{x)  —  J  J  d.r]dKA{'r),  cos(ija).  (4.1) 

Similarly,  the  complex  amplitude  perturbations  in  the  pupil  can  be  represented  by  the 
following  as  in  Eq.  (2.10): 

Pa(x)  ~  J  J  d.TjdKA(Ti,  K)e~^^^'^  sin(77Q').  (4.2) 

The  amplitude,  phase,  and  cross  correlations  were  developed  in  Chapter  2  by  com¬ 
puting  the  following 

Taa{^  =  {Pa{x)Pa(x-d!)),  (4.3) 

and 


Tpp{d)  =  {Ppix)Pp{S-^), 

and 


(4.4) 


Tapi^  =  {Pa{S)Ppix-^).  (4.5) 

where  the  angle  brackets  denote  the  expected  value  of  the  argument.  The  psds  were 
calculated  by  taking  the  Fourier  transforms  of  the  correlation  functions. 

The  next  step  will  be  to  generate  random  wavefront  screens  of  amplitude  and  phase 
with  the  proper  statistics.  The  screen  can  be  broken  into  a  phase  screen,  <^(x),  and  an 
amplitude  screen,  A{x).  The  two  can  be  combined  into  a  phase  and  amplitude  field  screen 


E{x)  =  (l  +  A{x))  (4.6) 

The  area  over  which  the  optical  field  perturbations  are  being  represented  by  the  FS 
is  a  square  of  dimension  Dp.  The  FS  representation  of  the  phase  and  amplitude  will  be 
periodic  in  both  the  x  and  y-directions  with  period  Dp.  The  FS  expansion  of  the  phase 
can  be  represented  as: 


OO  00 

^(®)=  E  E 

n=-oon'=-oo 


(4.7) 
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where  is  the  FS  coefficient  for  the  spatial  frequency  /  =  The  vectors  x 

and  y  are  unit  vectors  in  the  x  and  y-directions,  respectively,  and  x  and  y  are  the  x  and 
y  components  of  x.  The  coefficients  are  chosen  to  be  circular-complex  Gaussian, 
independent  random  variables  having  a  mean  square  value  given  by 


(4.8) 


where  E{-]  denotes  the  expected  value  of  the  argument  and  represents  the  phase  psd. 
We  require  Hermitian  symmetry,  =  c*-n,-n'-  require  the  zero  frequency  term 

to  be  a  real,  zero  mean  Gaussian  random  variable  [15].  In  Welsh’s  paper  he  used  the  geo¬ 
metric  phase  psd  in  the  above  equation,  but  one  could  just  as  easily  use  the  psds  developed 
for  the  phase  and  the  amplitude  pertmbations.  The  FS  expansion  of  the  amplitude  can  be 
represented  in  the  same  way  by  the  following: 


M^)=  *1  ^2^  (  ^  ^ 


n=— 00  n  — — 00 


Dr, 


n 

1). 


(4.9) 


where  c^J^|  is  the  FS  coefficient  for  the  spatial  frequency  /  =  ^  Similarly,  the 

coefficients  c^,^i  are  also  chosen  to  be  circular-complex  Gaussian,  independent  random 
variables  having  a  mean  square  value  given  by 


(4.10) 


and  due  to  Hermitian  symmetry,  Cn^n'  =  c*-n-ni-  The  zero  frequency  term  was  also  required 
to  be  a  real,  zero  mean  Gaussian  random  variable  [15]. 

We  could  simply  use  the  amplitude  and  phase  psds  to  calculate  independent  am¬ 
plitude  and  phase  effects,  but  since  they  are  not  independent  we  need  to  also  take  into 
account  their  cross  correlation.  Therefore,  we  need  to  satisfy  the  following  conditions: 


Dp’  Dp 


and 


and 


Dp’  Dp 


n  n 

Dp’  Dp 


(4.11) 


(4.12) 


(4.13) 
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To  achieve  these  conditions  we  must  satisfy  all  the  elements  in  the  correlation  matrix, 
i?c„  )  for  the  amplitude,  phase,  and  cross  correlations.  The  correlation  matrix  is  defined 
as 


(4.14) 


To  satisfy  Eqs.  (4.11),  (4.12),  and  (4.13)  consider 

(4.15) 

where  is  a  2  x  2  matrix  and  and  52  are  complex  zero  mean  unit  variance  independent 
Gaussian  random  numbers.  To  form  the  matrix  we  need  to  decompose  the  correla¬ 
tion  matrix,  J2c„  „(  j  so  that  the  phase,  amplitude  and  cross  correlation  are  all  satisfied. 
The  Cholesky  decomposition  is  performed  for  this  transformation,  since  i2c„  is  positive 
definite: 


.A 

‘'n.n' 


51 

52 


Rr  ,  =  C^C. 


Performing  the  decomposition  results  in  the  following  matrix 


CT  =  -^ 

D„ 


\j^AA  - 


Notice  what  happens  when  Rc^  is  computed: 


(4.16) 


(4.17) 


Rc  ,  =  E 

51 

[51  52]*  c 

52 
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c 


=  C^E 


91 

92 


52]* 


= 

= 


1  0 
0  1 


C  =  R 


(4.18) 


Therefore,  we  can  determine  what  random  numbers  to  use  for  the  screen  by  using 
Eq.  (4.15)  and  Eq.  (4.17): 


^n,n' 

1 

9l\/^<l»i> 

f.A 

^n,n' 

~  Dp 

(4.19) 


Notice  how  c^^j,  approaches  g^y/^AA  when  the  cross  correlation  approaches  zero  which  is 
the  desired  result.  Also,  the  phase  term  is  not  affected  by  the  cross  correlation. 

Welsh  approximates  the  phase  by  truncating  the  FS  expansion  in  Eq.  (4.7)  to  a  finite 
number  of  terms.  The  truncated  expansion  becomes  [15] 


^{S)  = 


N-l  N-1  (  / 

E  E 

-—(N-l)n'=-(N-l)  t  \  P  P^ 


=  2Re 


n=-{N-l)n'=-{N-l) 

TN-IN-1  (  /  /\ 


Ln=0  n'=0 


+  E  E  <n-»p|22-(^  +  ^)} 

n=l  „'=-(iV_l)  I  \^V  ^V)) 


(4.20) 


This  relationship  is  also  valid  for  the  amplitude,  A{x),  by  replacing  with 


4-2  Method  for  Screen  Generation 

Equation  (4.20)  can  be  expressed  in  matrix  algebra  notation,  allowing  for  simple 
implementation  on  mathematical  software  packages.  In  order  to  represent  the  equation  in 
matrix  algebra  form,  the  phase  screen  is  evaluated  on  a  M  x  M  square  grid  of  dimension  D, 
where  D  could  be  considered  the  diameter  of  the  optical  system  pupil.  In  most  instances 
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Dp  will  be  much  larger  than  D.  The  matrix  formulation  uses  the  following  vectors  in  the 
representation  of  the  phase  and  amplitude; 

X  :  column  vector  of  length  M  (x-coordinates  for  the  grid  points), 

y  :  column  vector  of  length  M  (y-coordinates  for  the  grid  points), 

f  :  column  vector  of  length  N  (spatial  frequency  with  nth  entry  being  [n  —  1]/I?p), 

1  :  column  vector  of  ones  of  length  M  [15]. 

Since  the  atmospheric  tmbulence  is  extended  over  the  entire  propagation  path  at 
constant  strength  of  turbulence,  the  amplitude  and  phase  screens  can  be  combined  and 
collapsed  into  one  screen  at  the  pupil.  The  phase  and  amplitude  can  be  represented  as  [15] 

X  +  y)  =  2Re  |[exp{j27r/a;  ^}] 

X  ((^^[exp  j27r/2/^]  +  ^^^[exp  -i27r/y^])  }  ,  (4.21) 


and 


A{x\'^ X  +  ly^ y)  =  2Re  |[exp{j27r/®  ^}] 

X  (C'i  [expi27r/i/^]  +  (^^[exp  -j27r/j/^])  }  ,  (4.22) 

where  Cf^  and  c*  are  N  X  N  matrices  containing  the  complex  random  variables  and 
and  are  N  x  N  matrices  containing  the  complex  random  variables  c^,^i . 

4.3  Field  Screen  Generation  Procedure 

The  psds  developed  in  Chapter  3  for  weah  turbulence  will  be  utilized  to  generate  the 
random  wavefront  phase  from  Eq.  (4.21)  and  the  random  amplitude  perturbation  from 
Eq.  (4.22).  The  following  considerations  and  steps  are  required  to  represent  the  phase  and 
amplitude  effectively. 

4.3.1  Set-up: 
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1.  FS  period,  D^,  Dp,  D^,  must  be  larger  than  the  non-negligible  support  of  Fpp,  Foa, 
and  Tap,  respectively.  If  Dp  is  not  selected  large  enough,  the  periodic  effects  of  Dp 
could  occur  in  the  structure  function  for  the  large  screen  case. 

2.  N  must  be  chosen  so  that  the  power  encompassed  in  the  simulated  field  sufiiciently 
represents  the  power  in  the  random  processes  ^{x)  and  A{x).  This  can  be  determined 
by  considering  the  ratio  of  the  power  in  ^{x)  to  the  power  in  ^(x)  and  the  ratio  of 
A{x)  to  the  power  in  A(x).  The  power  ratios  are  given  by 


E 

|<^(x)|2 

E 

mm 

for  the  phase. 


.n/d* 


(4.23) 


E 

E 

[|^(X)|2] 

J-N/D^ 


dp^AAi^ 


/_“  /_-  dp^AAi^ 


for  the  amplitude,  and 


(4.24) 


E 

A(x)^{x) 

E 

A{x)4){x) 

.N/D^  .N/D^ 

J-NinC 


-N/D'j, 


dp^A,j,{p) 


for  the  cross  term.  The  value  of  N/Dp  is  chosen  to  accurately  represent  the  wave¬ 
front’s  phase,  amplitude,  and  cross  statistics. 


3.  For  correlated  phase  and  amplitude  screens  Dp  as  well  as  N  must  be  the  same  for 
the  amplitude,  phase,  and  cross  statistics  (i.e..  Dp  =  D^  =  Dp  =  Dp)  to  properly 
represent  the  cross  statistics. 


4.  M  must  be  chosen  to  be  large  enough  to  sufficiently  sample  the  perturbations  in  the 
aperture  of  dimension,  D. 


4-3.2  Screen  Generation: 


1.  The  correlated  phase  and  amplitude  screens  will  be  calculated  by  filling  the  matrices 
C^,  C^,  C^,  and  with  circularly  complex  Gaussian  random  variables  where 
[s'!]*,*')  [ffaliji'!  and  are  complex  zero  mean  unit  variance  independent 
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Gaussian  random  numbers.  The  z,iTh  elements  of  the  matrices  will  contain  the 
following: 


(4.26) 


(4.27) 


Dp 


t— 1  i'— 1\ 
Dp  ’  Dp  y 

(i-\  i'-l\ 

V  Dp  >  Dp ; 

+  [53] 


^AA 


i  —  1  i'  —  1 

Dp  ’  Dp 


'^A4>  [Dp>  Dp  ) 
^ ) 


(4.28) 


^AA' 

f  i  —  1  1  — 

'i-l 

k  Dp  ’  Dp  ) 

' 

K  Dp  Dp  j 

G-1  l-i'’' 

^  Dp  ’  Dp  j 

(4.29) 


where  the  multiplicative  terms  represent  the  standard  deviation  of  the  random  num¬ 
bers,  [g2]i^,,  [53]*, and  [54]^^^,. 

2.  Replace  the  complex  random  values  of  the  first  row  and  column  of  and  with 
zeros. 

3.  Compute  the  random  screen,  E{x),  by  evaluating  Eq.  (4.21)  for  the  phase  and  Eq. 
(4.22)  for  the  amplitude  with  the  resulting  screen  having  the  following  form; 

E(x)=  (l  +  i(®))eJ^(^).  (4.30) 

4.  For  subsequent  independent  screen  realizations  go  to  step  1  under  subsection  4.3.2 
and  repeat  steps  1  through  3. 
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4-4  Numerical  Results 

In  this  section  the  specific  psds  calculated  in  Chapter  3  and  Eqs.  (4.21)  and  (4.22) 
will  be  used  to  generate  phase  and  amplitude  screen  realizations.  These  psds  are  valid  for 
a  plane  wave  using  a  constant  strength  of  turbulence  profile.  Monte  Carlo  experiments  for 
a  variety  of  conditions  are  performed  to  show  that  Eqs.  (4.21)  and  (4.22)  generate  screens 
with  the  correct  spatial  statistics. 


4-4-i  Weak  Turbulence  Atmospheric  Statistics.  Since  multiples  of  the  frequency 
1/Dp  are  used  in  the  FS  expansion  it  is  convenient  to  include  Dp  in  the  psds.  Therefore 
from  Eqs.  (3.16),  (3.17),  and  (3.18)  developed  in  the  previous  chapter  the  psds  become 


(  il^(27r) 

^aa(p)  =  167r®(0.033)(0.185)®/3-^^^^^^^ - - - 

(^^c| +(&)') 

['■  -  (2^) 


-11/3 

11/6 


=  0.01149- 


and 


^4><p{p)  =  IStt' 


'(0.033)(0.185)®/3-^^^^2J - P - 


X 


=  0.01149 


f'Dl 

\ro  j 

1  P 

2\  11/6 


1  + 


(4.32) 
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and 


( \  D^{iTx) 

=  16^^(0.033)(Q.185)V3  WW  ^  ^ 

('>^®l+(fe)) 

(tT^ 


wll/3  , 


sin^(7rp^AL) 


=  0.01149 


smr{T:p^\L) 


(4.33) 


4-4-^  Setup  for  independent  uncorrelated  phase  and  amplitude  screen  generation: 
Steps  1  and  2.  Prom  steps  1  and  2  outlined  in  subsection  4.3.1  the  FS  period,  Dp, 
and  the  number  of  FS  coefficients  N  must  be  selected  for  the  phase  and  amplitude.  For 
uncorrelated  screen  generation  the  cross  psd,  in  Eqs.  (4.28)  and  (4.29)  will  be  set 
to  zero.  Since  the  amplitude  and  phase  will  be  uncorrelated,  N  and  Dp  can  be  selected 
separately  to  best  represent  the  phase  and  amplitude  perturbations.  Prom  Fig.  2.3  it 
appears  that  a  value  of  D^  >  ALq  will  be  adequate  for  the  phase.  For  numerical  results 
shown  here  we  set  D^  —  ALq.  Using  Eq.  (4.23)  we  determine  that  N  >  15D^/Lo  is  required 
to  capture  99%  of  the  power  in  the  phase.  Since  we  chose  D^  =  ALg  at  least  iV  =  60  FS 
coefficients  are  required. 

Considering  the  amplitude,  from  Fig.  2.2  it  appears  that  a  value  of  Dp  >LolA  will 
be  adequate  for  the  amplitude.  For  numerical  results  shown  here  we  set  Dp  =  Lg/A.  Using 
Eq.  (4.24)  we  determine  that  N  >  720Dp/Lo  is  required  to  capture  99%  of  the  power  in 
the  amplitude.  Since  we  chose  Dp  =  Lg/A  at  least  N  =  180  FS  coefficients  are  required. 

We  also  need  to  consider  sample  spacing  of  the  screens.  In  these  numerical  results 
we  chose  N  =  ^D^/rg  for  the  phase  statistics.  In  the  following  results  we  chose  a  strength 
of  turbulence,  Lgjrg  =  100  which  for  the  phase  implies  that  N  =  (3/4)(£)^/io)(-^'oAo)  = 
1 . 4 . 100  =  300.  This  value  of  N  will  also  satisfy  the  power  requirements  for  the  amplitude 
using  the  Dp  =  LgfA. 

4-4-3  Simulation  results  for  independent  uncorrelated  amplitude  and  phase  screen 
generation.  Ten  thousand  random  screens  were  generated  in  a  Monte  Carlo  experiment 
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by  following  steps  1-4  of  subsection  4.3.2.  Each  screen  was  generated  on  a  square  grid  with 
a  screen  dimension  of  £>  =  lOOro  oi  D  =  \D^  for  the  phase.  The  large  screen  case  allowed 
us  to  compare  the  simulated  spatial  amplitude  and  phase  psds  to  the  theoretical  psds.  We 
chose  an  M  X  M  grid  for  the  screen  of  2D /to  x  2Dlro.  The  ensemble  based  wavefront 
phase  was  computed  to  be  187.61  square  radians,  which  compares  with  the  theoretical 
mean  square  value  (msv)  of  186.57  square  radians  and  the  amplitude  was  computed  to  be 
0.0877  which  compares  to  the  theoretical  msv  of  0.0879. 

Next,  a  screen  of  dimension  D  =  lOr^  is  considered.  This  time  one  thousand  screens 
were  generated  for  use  in  calculating  the  structure  functions.  An  M  x  M  grid  for  the 
field  screen  of  4D/r-o  X  AD /to  was  chosen.  Figure  4.1  shows  the  phase  structure  function 
calculated  for  the  ensemble  of  random  screens  along  with  the  theoretical  structure  function 
computed  from  Eq.  (2.38).  Figure  4.2  shows  the  calculated  amplitude  structure  function 
and  the  theoretical  function  computed  from  Eq.  (2.37).  The  ensemble  based  and  the 
theoretical  structure  functions  match  extremely  well.  Also,  notice  how  the  amplitude 
structure  function  away  from  the  center  of  the  screen  approaches  twice  the  mean  square 
value  of  the  amplitude  perturbations  calculated  in  the  previous  paragraph. 

Figiures  4.3  and  4.4  show  the  phase  and  amplitude  structure  functions  results  for  a 
larger  screen  dimension  of  D  =  40ro.  As  before  the  ensemble  based  and  the  theoretical 
structme  functions  match  extremely  well. 

Finally,  the  amplitude  and  phase  structure  functions  are  calculated  for  a  very  large 
screen  dimension  of  I?  =  lOOrQ.  For  this  case,  Dp/Lo  for  the  amplitude  must  be  selected 
as  Dp /Lo  >  5/8  in  order  to  keep  the  periodic  nature  of  Dp  from  corrupting  the  structure 
function.  The  phase  structure  function  is  in  Fig.  4.5  and  matches  very  well  with  the 
theoretical  function.  If  N  =  300  is  selected,  the  amplitude  psd  will  be  under-represented, 
since  N  >  720Dp/Lo  will  no  longer  be  satisified.  Therefore,  if  we  let  N  =  450,  over  99  % 
of  the  power  in  the  amplitude  psd  will  be  captured.  Figure  4.6  shows  that  the  resulting 
amplitude  structure  function  is  slightly  under  represented  but  within  1%  of  the  theoretical 
plot. 
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Figure  4.1  Theoretical  and  simulated  phase  structure  function,  Dpp{S)  plotted  versus 
|d|/Lo  for  a  von  Karman  turbulence  power  spectrum  with  D^/Lq  =  4, 
D^/tq  =  400,  L/Lo  =  1000,  and  XjLo  =  5  x  10“®.  The  screen  dimension  is 
D  =  lOro  and  the  number  of  FS  coefficients  is  N  =  300. 
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Figure  4.2  Theoretical  and  simulated  amplitude  structure  function,  Daa{^  plotted  ver¬ 
sus  |dj/Lo  for  a  von  Karman  tiurbulence  power  spectrum  with  D^/Lo  =  1/4, 
Dpiro  =  25,  LjLo  =  1000,  and  A/Lp  =  5  x  10~®.  The  screen  dimension  is 
D  =  lOro  and  the  number  of  FS  coefficients  is  iV  =  300. 
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Figure  4.3  Theoretical  and  simulated  phase  structure  function,  Dpp{S),  plotted  ver¬ 
sus  |d|/io  for  a  von  Karman  tiurbulence  power  spectrum  with  D^/Lg  =  4, 
D^/tq  =  400,  LjLo  =  1000,  and  XJLo  =  5  x  10“®.  The  screen  dimension  is 
D  =  40ro  and  the  number  of  FS  coefficients  is  =  300. 


-0.2  -0.16  -0.1  -0.05  0  0.05  0.1  0.16  0.2 


Figure  4.4  Theoretical  and  simulated  amplitude  structure  function,  Daa{^,  plotted  ver¬ 
sus  \d\/Lo  for  a  von  Karman  turbulence  power  spectrum  with  Dp /Lg  =  1/4, 
Dp/vg  =  25,  L/Lg  =  1000,  and  X/Lg  =  5  x  10“®.  The  screen  dimension  is 
D  =  40ro  and  the  number  of  FS  coefficients  is  N  =  300. 
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Figure  4.5  Theoretical  and  simulated  phase  structure  function,  Dpp{S),  plotted  ver¬ 
sus  \d\lLo  for  a  von  Karman  turbulence  power  spectrum  with  D^/Lq  =  4, 
D^/ro  =  400,  LjLo  =  1000,  and  X/Lq  =  5  x  10~®.  The  screen  dimension  is 
D  =  lOOfo  and  the  number  of  FS  coefficients  is  =  300. 
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Theoretical  and  simulated  amplitude  structure  function,  Daai^,  plotted  ver¬ 
sus  \d\/Lo  for  a  von  Karman  turbulence  power  spectrum  with  Dp/Lg  =  5/8, 
Dplvo  =  62.5,  LjLo  =  1000,  and  X/Lg  =  5  x  10~®.  The  screen  dimension  is 
D  =  lOOro  and  the  number  of  FS  coefficients  is  AT  =  450. 
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4-4-4  Setup  for  independent  correlated  phase  and  amplitude  screen  generation: 
Steps  1,  2,  and  3.  Prom  steps  1  and  2  outlined  in  subsection  4.3.1  the  FS  period,  Dp, 
and  the  number  of  FS  coeflScients  N  must  be  selected  for  the  phase,  amplitude,  and  cross 
psds.  Correlated  phase  and  amplitude  screen  generation  requires  the  same  values  for  Dp 
and  N.  Prom  subsection  4.4.2  it  is  apparent  that  the  amplitude  and  phase  requirements 
for  Dp  and  N  are  not  the  same.  First  Dp  must  be  selected,  since  it  must  be  chosen  large 
enough  to  represent  the  amplitude,  phase,  and  cross  statistics  sufficiently.  In  considering 
the  cross  statistics,  from  Fig.  2.4  it  appears  that  a  value  of  Dp  >  Lq  will  be  adequate 
for  the  cross  correlation.  Prom  Eq.  (4.25)  N  >  lOODpfLo  is  determined  to  be  required 
to  capture  99%  of  the  power  in  the  cross  statistics.  In  order  to  sufficiently  represent  the 
phase,  amplitude,  and  cross  statistics  Dp  =  1.5  is  chosen.  To  captiure  99%  of  the  power  in 
the  amplitude,  N  >  720Dp/Lo  must  be  satisfied,  requiring  N  >  1080.  This  requirement 
is  much  higher  than  the  practical  limits  of  the  computation  speed,  especially  given  the 
number  of  iterations  in  the  Monte  Carlo  experiments.  N  =  400  FS  coefficients  was  chosen 
which  resulted  in  a  representation  of  over  99%  in  the  phase  and  cross  statistics  and  92.3% 
for  the  amplitude. 

4.4-5  Simulation  results  for  independent  correlated  amplitude  and  phase  screen 
generation.  Ten  thousand  random  screens  were  generated  in  a  Monte  Carlo  experiment 
by  following  steps  1-4  of  subsection  4.3.2.  Each  screen  was  generated  on  a  square  grid  with  a 
screen  dimension  of  Z?  =  lOOro.  The  large  screen  case  allowed  us  to  compare  the  simulated 
spatial  amplitude,  phase,  and  cross  psds  to  the  theoretical  psds.  As  before,  we  chose  an 
MxM  grid  for  the  screen  of  2£)/roX2jD/ro.  Using  £>p  =  1.5  and  iV  =  400  determined  in  the 
previous  section  the  ensemble  based  phase,  amplitude  and  cross  values  were  determined. 
The  ensemble  based  wavefront  phase  was  computed  to  be  185.41  square  radians,  which 
compares  with  the  theoretical  mean  square  value  (msv)  of  186.57  square  radians  and  the 
amplitude  was  computed  to  be  0.0812  which  compares  to  the  theoretical  msv  of  0.0879.  The 
ensemble  based  mean  cross  value  of  the  amplitude  and  phase  perturbations  was  computed 
to  be  0.2416  which  compares  with  the  theoretical  value  of  0.2418. 
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Next  a  screen  of  dimension  D  =  lOrg  is  considered.  Ten  thousand  screens  are 
generated  to  calculate  the  structure  functions.  Again,  an  M  x  M  grid  for  the  field  of 
AD/to  X  AD/vo  is  selected.  Figure  4.7  shows  the  phase  structure  function  calculated  for 
the  ensemble  of  random  screens  along  with  the  theoretical  structure  function  computed 
from  Eq.  (2.38).  Figure  4.8  shows  the  calculated  amplitude  structure  function  and  the 
theoretical  function  computed  from  Eq.  (2.37).  Figure  4.9  shows  the  calculated  cross 
structure  function  and  the  theoretical  function  computed  from  Eq.  (2.39).  The  ensem¬ 
ble  based  and  the  theoretical  structure  functions  match  extremely  well  for  the  phase  and 
cross,  but,  as  expected,  the  simulated  amplitude  structure  function  drops  slightly  below 
the  theoretical  plot. 


Figure  4.7 


Theoretical  and  simulated  phase  structure  function,  Dpp{S),  plotted  versus 
|d|/Xo  for  a  von  Karman  turbulence  power  spectrum  with  Dp/Lo  =  1.5, 
Dp/ro  =  150,  L/Lo  =  1000,  and  X/Lo  =  5  x  10“®.  The  screen  dimension  is 
D  =  lOr-o  and  the  number  of  FS  coefficients  is  N  =  400.  The  phase  screens 
used  in  this  calculation  were  correlated  with  the  amplitude  screens  by  the 
cross  correlation. 


Figures  4.10,  4.11,  4.12  show  the  phase,  amplitude,  and  cross  structure  functions 
results  for  the  very  large  screen  dimension  of  D  =  IOOt-q.  As  before,  the  ensemble  based 
and  the  theoretical  structure  functions  match  extremely  well  for  the  phase  and  cross,  but 
the  ensemble  based  amplitude  structure  function  drops  slightly  below  the  theoretical. 
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Figure  4.8 
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Theoretical  and  simulated  amplitude  structure  function,  Daai^,  plotted  ver¬ 
sus  \d\/Lo  for  a  von  Karman  turbulence  power  spectrum  with  Dp/Lo  =  1.5, 
Dp/ro  =  150,  LjLo  =  1000,  and  A/To  =  5  X  10“®.  The  screen  dimension 
is  D  =  IOt-q  and  the  number  of  FS  coefficients  is  iV  =  400.  The  amplitude 
screens  used  in  this  calculation  were  correlated  with  the  phase  screens  by  the 
cross  correlation. 


Figure  4.9  Theoretical  and  simulated  cross  structure  function,  Dap{^,  plotted  versus 
|d|/To  for  a  von  Karman  turbulence  power  spectrum  with  Dp/Lo  =  1.5, 
Dp/vo  =  150,  L/Lo  =  1000,  and  A/Tq  =  5  x  10“®.  The  screen  dimension  is 
D  =  lO^o  and  the  number  of  FS  coefficients  is  N  =  400. 
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Figure  4.10 


Figure  4.11 


Theoretical  and  simulated  phase  structure  function,  Dpp{S),  plotted  versus 
\d\lLo  for  a  von  Karman  turbulence  power  spectrum  with  Dp/Lo  =  1.5, 
Dp/ro  =  150,  LjLo  =  1000,  and  X/Lq  =  5  x  10“®.  The  screen  dimension  is 
D  =  lOOro  and  the  number  of  FS  coefficients  is  =  400.  The  phase  screens 
used  in  this  calculation  were  correlated  with  the  amplitude  screens  by  the 
cross  correlation. 
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Theoretical  and  simulated  amplitude  structure  function,  Daa{^,  plotted 
versus  \d\/Lo  for  a  von  Karman  tmrbulence  power  spectrum  with  Dp/Lo  = 
1.5,  Dp/ro  =  150,  L/Lo  =  1000,  and  X/Lq  =  5  x  10~®.  The  screen  dimension 
is  D  =  lOOr-Q  and  the  number  of  FS  coefficients  is  N  =  400.  The  amplitude 
screens  used  in  this  calculation  were  correlated  with  the  phase  screens  by 
the  cross  correlation. 
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Figure  4.12 


Theoretical  and  simulated  cross  structure  function,  Dap{^,  plotted  versus 
\d\/Lo  for  a  von  Karman  turbulence  power  spectrum  with  Dp/Lo  =  1.5, 
Dp/vo  =  150,  LjLo  =  1000,  and  XJLq  =  5  x  10“®.  The  screen  dimension  is 
D  =  lOOro  and  the  number  of  FS  coefficients  is  iV  =  400. 
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Figure  4.13  shows  a  random  phase  screen  of  dimension  D  =  lOr-p  and  Fig.  4.14 
shows  the  correlated  random  amplitude  screen.  Notice  how  much  larger  the  phase  screen 
magnitude  is  than  the  amplitude  screen.  Figs.  2.5  and  2.6  from  Chapter  2  illustrated  the 
mean  square  value  of  the  amplitude  perturbations  will  get  larger  and  the  phase  will  get 
smaller  as  decreases  for  constant  Lo/to.  Therefore,  if  L/Lg  or  A/Tq  were  larger  the 
magnitude  of  the  amplitude  perturbations  would  be  higher  and  the  magnitude  of  the  phase 
perturbations  would  be  lower.  It  is  also  evident  that  the  amplitude  perturbations  have  a 
higher  spatial  frequency  content  than  the  phase  perturbations  which  can  be  predicted  from 
the  power  spectral  density  plot  in  Fig.  3.1. 


Figure  4.13  Random  phase  screen,  ^(d),  plotted  versus  \d\lLo  for  a  von  Karman  tur¬ 
bulence  power  spectrum  with  DpjLo  =  4,  Dplvg  =  400,  LjLo  =  1000,  and 
X/Lo  =  5  X  10“®.  The  screen  dimension  is  D  =  lOr,,  and  the  number  of  FS 
coefficients  is  N  =  300. 
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Figure  4.14 


Random  amplitude  screen,  plotted  versus  |d|/Lo  for  a  von  Karman 

turbulence  power  spectrum  with  Dp/Lp  =  1/4,  Dpjro  =  25,  L/Lo  =  1000, 
and  X/Lo  =  5  x  10“®.  The  screen  dimension  is  I?  =  lOr,,  and  the  number 
of  FS  coefiicients  is  iV  =  300. 
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V.  Amplitude  Effects  on  Full  Aperture  Tilt  Sensing 

The  random  tilt  over  the  aperture  of  an  imaging  system  has  the  effect  of  moving  an  image 
around  in  the  image  plane  of  an  optical  system.  Long  exposure  images  will  effectively 
average  this  movement  over  a  long  period  of  time  blurring  the  image  and  reducing  the 
resolution  similar  to  a  convolution  operation.  A  short  exposure  image  essentially  freezes 
the  random  tilt  allowing  for  much  higher  resolution.  The  tip-tilt  mirror  corrects  the  overall 
mean  tilt  over  the  aperture  of  the  wavefront,  significantly  reducing  the  movement  of  the 
image  in  the  image  plane.  It  is  simply  a  fiat  mirror  that  can  be  adjusted  in  two  orthogonal 
directions.  Adaptive  optics  systems  began  with  a  tilt  correction  system  since  the  random 
tilt  contains  87%  of  the  power  of  the  wavefront  phase  fluctuations.  In  modern  adaptive 
optics  systems  the  tip-tilt  mirror  stage  separately  removes  the  random  tilt  in  a  wavefront 
reducing  the  dynamic  range  requirements  in  the  fiexible  or  faceted  deformable  mirror. 
Refer  back  to  Fig.  1.1  for  a  simplified  adaptive  optics  system  schematic.  Dynamic  range 
refers  to  the  deflection  magnitude  of  each  actuator  or  facet  of  the  mirror  required  to  correct 
the  higher  order  phase  errors  [12]. 

5.1  Full  aperture  tilt  estimation  using  image  centroid 

There  are  a  number  of  algorithms  and  systems  used  to  measure  the  overall  tilt  of 
the  wavefront  in  the  aperture.  One  such  measurement  tracks  the  movement  of  the  image 
centroid  in  the  focal  plane  of  a  full  aperture  lens.  Figure  5.1  illustrates  the  fundamental 
structure  of  the  centroid  tracker. 

When  the  centroid  of  the  image  is  located  in  the  center  of  the  focal  plane  the  overall 
tilt  of  the  phase  perturbations  in  the  aperture  is  zero.  Using  the  Fourier  transforming 
property  of  a  lens  the  image  intensity,  I{u,v),  can  be  determined  from  the  optical  field, 
E{x,y),  by  the  following: 

=y^{El^,y)}r{Elx,y)r,  (5.1) 

where  T  is  the  Fourier  transform  operator,  A  is  the  average  wavelength  of  the  light,  /  is 
the  focal  length  of  the  imaging  lens,  (a:,y)  are  the  coordinates  in  the  object  plane,  and 
are  the  coordinates  in  the  focal  plane  [7]. 
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Figure  5.1  Simplified  configuration  of  a  full  aperture  tilt  centroid  tracker.  {Cx,Cy)  is 
the  location  of  the  image  centroid  and  /  is  the  focal  length  of  the  lens. 

The  centroid  of  the  image  can  be  determined  by  performing  a  discrete  first  moment 
calculation  to  determine  the  x  and  y  locations  of  the  centroid: 

k 

4  =  ^4 -  (5.2) 

i=l 

and 

k 

4  ~  ’  (^•^) 

i—1 

where  Mj  is  the  relative  irradiance  on  each  element  of  the  detector  in  the  focal  plane  and 
{xi,yi)  are  the  locations  of  each  element.  {Ix,Iy)  represents  the  discrete  location  of  the 
image  centroid  on  the  detector. 

Simulations  will  be  run  to  determine  the  image  in  the  focal  plane  by  evaluating  Eq. 
(5.1)  and  using  the  FFT  to  compute  the  Fourier  transform  of  the  wavefront.  The  centroid 
of  the  image  will  be  determined  and  related  to  the  slope  of  the  wavefront  in  the  aperture 
of  the  system  from  the  relationships  in  Eqs.  (5.1),  (5.2),  and  (5.3).  Since  the  FFT  will  be 
implemented,  a  relation  between  the  sample  size  in  the  apertme  to  the  sample  size  of  the 
focal  plane  detector  is  required.  For  a  m  X  m  matrix  used  in  the  computation  of  the  image 
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the  following  relationships  relate  the  Foirrier  frequency  sample  size  to  the  sample  size  in 
the  space  domain: 


^fx  = 


mAx 


=  7;:Ay 

where  Afx  and  Afy  are  the  sample  size  in  the  Fourier  domain  and  Ax  and  Ay  are  the 
sample  size  in  the  space  domain  [5].  Now,  using  Eq.  (5.1)  one  can  determine  the  sample 
size  in  the  image  plane 


Au  = 


mAx 


Av  = 


mAy^ 


where  Au  and  Av  are  sample  size  in  the  focal  or  detector  plane.  Prom  Eqs.  (5.2)  and 
(5.3)  one  can  determine  the  actual  location  of  the  centroid  by  computing  C®  =  IxAu  and 

Cy  =  IxAV, 

Now,  consider  how  the  centroid  location  in  the  focal  plane  is  related  to  the  wavefront 
slope  over  the  aperture.  Prom  Fig.  5.1  the  location  of  the  image  centroid,  {Cx,  Cy),  is 
related  to  the  orthogonal  components  of  the  wavefront  slope  s  by 

_ _  ZTrCa 

St.  —  ^  -  (5.8) 


_  7  Cy  _  2TrCy  /r  n\ 

~  ~  !  (^•^) 

where  (sx,  Sy)  represents  the  x-  and  y-directed  components  of  the  vector  s  [12].  Expressing 
the  discrete  representation  of  the  centroid  the  wavefront  slope  can  be  simplified  by 


2'kIxAu  _  2'kIx  a/ 
A/  A/  mAx 
27r4 
mAx 


(5.10) 
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and 


_  2ttIxAv  _  2'Kly  A/ 

A/  A/  mAy 
27r4 
mAy' 

Recall  that  m  x  m  is  the  dimension  of  the  matrix  in  the  FFT  calculation. 


(5.11) 


5.1.1  Tilt  estimation  calculation  using  random  screens.  The  random  phase  and 
amplitude  screen  generator  outlined  in  Chapter  4  will  be  used  to  determine  the  effect 
the  amplitude  perturbations  have  on  the  wavefront  slope  estimation,  (si,  Sy).  The  random 
screen  consists  of  a  M  X  M  matrix  of  dimension,  D.  Therefore,  Ax  =  D/M  and  Ay  =  D/M 
and  (sj;,  Sy)  becomes 


2TrIx 

m.D 


(5.12) 


and 


2'Kly 


(5.13) 


Using  Eqs.  (5.6)  and  (5.7)  one  can  determine  the  size  of  the  FFT  matrix  that  must 
be  used  to  determine  the  image  in  the  focal  plane  given  the  particular  sample  size  of  the 
detector,  Au  and  the  optical  parameters  A,  /,  and  D.  Solving  for  m,  Eq.  (5.6)  becomes 


_  A/  _  A/M  _  A/#M 
^  AxAu  DAu  Au  ’ 

where  Au  =  Av,  Ax  =  Ay,  and  /#  =  f/D  is  the  f-number  of  the  lens. 


(5.14) 


5.2  Simulations  using  the  random  screen  generator  with  uncorrelated  amplitude  and 
phase  screens 

The  screen  generator  will  be  used  to  investigate  how  the  amplitude  perturbations 
affect  the  tilt  calculation.  The  mean  square  error  in  the  tilt  calculation  will  be  determined 
for  an  ensemble  of  square  screens.  The  estimated  tilt  will  be  subtracted  from  the  phase 
screens  to  determine  how  well  the  algorithm  performs  for  different  situations. 
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Turbulence  and  Optical  Parameters 

Parameter 

Range 

Default 

low  high 

Propagation  length  over  outer  scale,  L/Lq. 

1000-50,000 

10,000 

Strength  of  turbulence,  Lo/rg. 

50-500 

200 

Wavelength  over  outer  scale,  X/Lg. 

Diameter  over  Fried  parameter,  D/rg. 

2-100 

10 

Table  5.1  Atmospheric  turbulence  and  optical  parameters  used  in  simulations. 


The  amplitude  and  phase  screens  used  in  these  simulations  are  uncorrelated.  Section 
5.3  provides  justification  for  using  the  uncorrelated  screens  over  the  correlated  screens. 

The  following  parameters  Djro,  LfLo,  LqIto.,  and  XJLq  as  shown  in  Table  5.1  will 
be  studied  to  determine  their  effect  on  the  tilt  estimation.  Square  screens  of  width,  D,  will 
be  used  to  simulate  atmospheric  tmbulence  over  a  square  aperture.  As  each  parameter  is 
varied,  all  other  parameters  will  remain  constant.  Since  many  of  the  parameters  are  inter¬ 
related,  care  must  be  taken  in  the  analysis  of  the  results.  For  instance,  normally  decreasing 
LjLo  will  result  in  less  turbulence  effects.  Since  Dfro  is  proportional  to  the  integrated 
turbulence  profile  and  held  constant,  decreasing  LjLo  had  the  effect  of  simulating  shorter 
but  more  concentrated  turbulence  strength  profiles. 


5.2.1  Mean  square  error  in  the  tilt  estimation.  The  estimated  tilt  will  be  cal¬ 
culated  for  each  random  screen,  the  square  error  in  the  estimation  determined  for  each 
realization,  and  the  mean  calculated  for  the  entire  ensemble.  The  mean  square  error  in  the 
tilt  estimation  can  be  defined  as 


Tmse 


{{tx  +  {ty  Sy)^) 

2 


(5.15) 


where  tx  and  ty  are  the  x—  and  y— directed  components  of  the  actual  tilt  and  the  angle 
brackets  represent  the  ensemble  average  operation.  The  Tmse  for  a  random  amplitude  and 
phase  screen  will  be  compared  with  the  error  in  the  tilt  estimation  using  a  phase  only 
screen.  The  components  of  the  actual  tilt,  tx  and  ty,  are  calculated  by  taking  the  average 
slope  over  the  phase  screen.  The  mean  square  tilt  of  the  phase  screens  is  defined  as 


(5.16) 
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The  standard  deviation  of  the  mean  for  the  Monte  Carlo  experiments  can  be  deter¬ 


mined  by 


Vn~  Vn’ 


(5.17) 


where  a  is  the  actual  standard  deviation  of  the  population,  s  is  the  standard  deviation  of 
the  sample,  and  N  is  the  number  of  iterations  in  the  simulation  [1].  Using  Eq.  (5.17), 
an  experiment  of  1200  random  screens  or  iterations  was  determined  to  be  sufficient  in  the 
Monte  Carlo  calculations  of  the  tilt.  The  standard  deviation  in  the  mean  for  Tmse  was 
approximately  l/35th  of  the  mean. 


In  the  Monte  Carlo  experiments  the  following  calculations  with  performed  for  each 
iteration.  First,  the  actual  tilt  is  calculated  from  the  phase  screen.  Then,  the  full  aperture 
tilt  estimation  algorithm  is  performed  using  Eqs.  (5.12)  and  (5.13)  on  the  optical  field 
screen  consisting  of  the  phase  and  amplitude  perturbations.  Next,  we  use  Eqs.  (5.12)  and 
(5.13)  again,  but  this  time  we  include  only  the  phase  perturbations  in  the  screen.  Finally, 
the  mean  square  error  in  the  tilt,  Tmaei  is  calculated  for  the  phase  only  screen  and  the  phase 
and  amplitude  screen  using  Eq.  (5.15).  The  difference  between  the  represents  the 

effect  the  amplitude  perturbations  have  on  the  tilt  estimation.  The  mean  square  tilt,  Tms, 
and  Tmse  ar®  presented  in  the  plots  in  radians^  per  number  of  Tq  squared  (or  radVg^). 
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In  the  first  case,  the  propagation  length  over  outer  scale,  L/Lq,  was  varied.  The 
mean  square  error  in  the  tilt  estimation,  T^sej  is  plotted  in  Fig.  5.2  for  both  cases. 
Normally  as  the  propagation  path  lengthens  the  turbulence  effects  become  more  severe, 
but  since  Djro  and  Lo/to  are  held  constant,  only  the  proportion  of  the  amplitude  to 
the  phase  perturbations  changes.  Figs.  2.5  and  2.6  from  Chapter  2  show  that  the  mean 
square  value  of  the  amplitude  increases  with  L/Lg  and  the  mean  square  value  of  the  phase 
decreases  when  all  other  factors  are  held  constant.  Therefore  as  LjLo  gets  larger,  the 
larger  amplitude  perturbations  cause  the  Tmae  for  the  phase  and  amplitude  screen  case 
to  increase.  Whereas,  the  phase  only  Tmse  decreases  as  LjLo  increases  due  to  the  slight 
decrease  in  phase  perturbations.  In  addition,  as  these  phase  perturbations  decrease  with 
smaller  L/Lg  the  mean  square  tilt,  T^s,  decreases  as  expected. 


Figure  5.2  The  mean  square  error  in  the  tilt  estimation,  Tmsei  and  the  mean  square 
tilt,  Tms,  are  plotted  versus  LjLg  for  D/vg  =  10,  Lg/rg  =  200,  and  X/Lg  = 
5.0x10-8.  The  Tmse  is  calculated  twice,  once  using  the  phase  and  amplitude 
screen  case  and  once  for  the  phase  screen  only  case.  The  Tms  is  calculated 
by  computing  the  average  slope  of  the  phase  over  the  screen. 
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Next  the  strength  of  turbulence,  Lo/ro,  was  varied  and  the  Tms  and  Tmse  were  calcu¬ 
lated  again.  The  mean  square  error  for  both  cases  and  the  mean  square  tilt  are  shown  in 
Fig.  5.3.  As  LoItq  increases  and  D/to  is  held  constant  the  lower  order  phase  perturbations 
(e.g.,  tilt)  increase  while  the  higher  order  perturbations  decrease.  The  upward  trend  in 
the  tilt,  Tmsi  is  due  to  larger  outer  scales,  X-o,  contributing  to  higher  tilt  in  the  wavefront, 
but  flatter  smoother  variations.  The  mean  square  value  of  the  amplitude  perturbations 
increase  with  the  strength  of  turbulence  as  shown  in  Fig.  2.6.  This  increase  in  ampli¬ 
tude  effects  causes  greater  error  in  the  tilt  estimate  for  the  phase  and  amplitude  screen. 
Whereas,  Tmae  for  the  phase  only  case  decreases  as  Lo/to  increases,  since  the  tilt  estimator 
performs  better  for  flatter  phase  screens. 


To 


Figure  5.3  The  mean  square  error  in  the  tilt  estimation,  Trasei  and  the  mean  square  tilt, 
Tmsi  are  plotted  versus  Lo/to  for  Djro  =  10,  L/Lq  =  10000,  and  X/Lq  = 
5.0x10-8.  The  Tfnse  is  calculated  twice,  once  using  the  phase  and  amplitude 
screen  case  and  once  for  the  phase  screen  only  case.  The  Tms  is  calculated 
by  computing  the  average  slope  of  the  phase  over  the  screen. 
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Next,  the  tilt  estimation  is  investigated  for  changes  in  XjLo  in  Fig.  5.4.  Similar 
to  LjLoi  Figs.  2.5  and  2.6  from  Chapter  2  also  show  that  the  mean  square  value  of  the 
amplitude  increases  with  A/Lq  and  the  mean  square  value  of  the  phase  decreases  when 
the  other  factors  are  held  constant.  As  XfLo  gets  larger,  Tmse  for  the  phase/amplitude 
case  increases  due  to  the  larger  amplitude  perturbations.  Whereas,  the  phase  only  Tmse 
decreases  slightly  as  X/Lo  increases  due  to  the  decrease  in  phase  pertmbations.  In  addition, 
as  these  phase  perturbations  decrease  the  mean  square  tilt,  T^s,  decreases  slightly  as 
expected. 


Figure  5.4  The  mean  square  error  in  the  tilt  estimation,  T^sej  and  the  mean  square  tilt, 
Tms,  are  plotted  versus  XfLo  for  D/ro  =  10,  LfLo  =  10000,  and  Lo/to  =  200. 
The  Tjnse  is  calculated  twice,  once  using  the  phase  and  amplitude  screen  case 
and  once  for  the  phase  screen  only  case.  The  Tma  is  calculated  by  computing 
the  average  slope  of  the  phase  over  the  screen. 
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Finally,  the  tilt  estimation  is  investigated  for  changes  in  D/to  in  Fig.  5.5.  Since 
L/Lo,  Lo/ro,  and  A/Lq  are  held  constant  the  strength  and  nature  of  the  phase  and  ampli¬ 
tude  perturbations  remain  relatively  constant  over  a  fixed  area.  Over  very  large  apertures 
the  phase  tilt  is  very  small,  but  as  the  aperture  size  decreases  the  mean  square  tilt,  Tms, 
increases.  Therefore,  as  D/ro  gets  larger,  the  tilt  over  the  entire  aperture  decreases  signif¬ 
icantly.  The  error,  Tmsei  decreases  for  both  cases  as  Tms  decreases. 


Figure  5.5  The  mean  square  error  in  the  tilt  estimation,  Tmse^  and  the  mean  square 
tilt,  Tmsi  are  plotted  versus  D/to  for  A/Tp  =  5.0  X  10“®,  L/Lq  =  10000, 
and  Lot  To  =  200.  The  Tmse  is  calculated  twice,  once  using  the  phase  and 
amplitude  screen  case  and  once  for  the  phase  screen  only  case.  The  Tms  is 
calculated  by  computing  the  average  slope  of  the  phase  over  the  screen. 
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5.2.2  Strehl  ratio  calculation.  The  Strehl  ratio  will  be  calculated  before  and 
after  removing  the  tilt  in  the  screen  for  the  ensemble  of  random  screens.  The  ensemble  of 
screens  will  be  generated  in  a  Monte  Carlo  experiment.  The  point  spread  function  is  the 
image-plane  intensity  distribution  resulting  from  imaging  a  point  source  [12].  The  point 
spread  function  (psf)  can  be  represented  as 

s{u,v)  =  *  ^(0, 0) 

=  I{u,v),  (5.18) 

where  6(0, 0)  is  the  Dirac  delta  function  and  *  is  the  convolution  operator.  The  psf  will 
be  calculated  for  each  screen  and  averaged  over  the  ensemble  to  simulate  a  long  exposure 
image.  Each  corrected  screen  will  be  determined  by  subtracting  the  estimated  tilt  from 
the  phase  screen.  The  Strehl  ratio  is  defined  zis  the  ratio  of  maximums  between  the  actual 
psf  and  the  diffraction  limited  psf.  The  Strehl  ratio  will  be  calculated  for  fom:  different 
cases: 


1.  The  actual  tilt,  tx  and  ty,  will  be  determined  by  taking  the  average  slope  over  the 
entire  phase  screen.  The  actual  tilt  will  be  subtracted  from  the  phase  screen  and 
the  Strehl  ratio,  SRut,  will  be  calculated  for  the  ensemble  average  psf  of  the  cor¬ 
rected  phase  and  original  amplitude  screens.  The  subscript,  UT,  refers  to  the  screens 
“untilted”  by  the  actual  tilt. 

2.  The  estimated  tilt,  Sx  and  Sy,  will  be  determined  from  Eqs.  (5.1),  (5.12),  and  (5.13) 
by  using  the  centroid  of  the  phase  screen  only  image.  The  estimated  tilt  will  be 
subtracted  from  the  phase  screen  and  the  Strehl  ratio,  SRutPi  will  be  calculated  for 
the  ensemble  average  psf  of  the  corrected  phase  and  the  original  amplitude  screens. 
The  subscript,  UTP,  refers  to  the  screens  “untilted”  by  the  estimated  tilt  from  the 
phase  screen  image  centroid. 

3.  The  estimated  tilt,  Sx  and  Sy,  will  be  determined  from  Eqs.  (5.1),  (5.12),  and  (5.13) 
by  using  the  centroid  of  the  phase  and  amplitude  screen  image.  The  estimated 
tilt  will  be  subtracted  from  the  phase  screen  and  the  Strehl  ratio,  SRutaPi  will 
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be  calculated  for  the  ensemble  average  psf  of  the  corrected  phase  and  the  original 
amplitude  screens.  The  subscript,  UTAP,  refers  to  the  screens  “untilted”  by  the 
estimated  tilt  from  the  amplitude/phase  screen  image  centroid. 

4.  The  Strehl  ratio,  SR,  will  be  calculated  for  the  uncorrected  phase  screen  using  the 
ensemble  average  psf  of  the  combined  phase  and  amplitude  screens. 

Strehl  ratio  calculations  1  through  3  will  have  the  effect  of  decreasing  the  movement  of 
the  image  in  the  focal  plane.  Whereas,  the  uncorrected  phase  screen  case  of  calculation 
4  above  will  be  averaged  over  a  wider  range  of  image  motion.  The  corrected  screen  with 
the  least  amount  of  image  motion  over  the  ensemble  will  result  in  the  highest  Strehl  ratio. 
Figure  5.6  shows  a  block  diagram  of  the  fom:  different  Strehl  ratio  calculations  performed. 
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Random  screen  with  phase 
perturbations  oniy 


Random  screen  with  phase  & 
ampiitude  perturbations 


Figure  5.6  Block  diagram  of  Strehl  ratio  calculation  procedure.  Random  screens  are 
generated  by  following  the  Screen  Generation  procedure  in  Sec.  4.3.2.  Strehl 
ratios  SRut,  SRutp,  and  SRutap  are  for  screens  compensated  by  the  actual 
tilt,  estimated  tilt  using  phase  screen  image  centroid,  and  estimated  tilt  using 
phase/amplitude  screen  image  centroid,  respectively.  The  fourth  Strehl  ratio, 
SR,  is  calculated  for  the  random  uncompensated  phase/amplitude  screen. 
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Figure  5.7  shows  the  point  spread  function  of  the  diffraction  limited  unaberated 
screen.  The  uncorrected  random  screen’s  psf  for  the  short  exposure  image  is  shown  in 
Fig.  5.8.  The  simulations  were  run  by  performing  a  Monte  Carlo  experiment  to  generate 
random  screens.  By  taking  the  ensemble  mean  of  the  psf ’s  of  each  screen  the  long  exposure 
point  spread  functions  were  determined.  The  uncorrected  long  exposure  psf  is  shown  in 
Fig.  5.9  (case  4  in  Fig.  5.6).  The  long  exposure  psf  for  the  screens  compensated  by  the 
phase  only  image  centroid  estimated  tilt  is  shown  in  Fig.  5.10  (case  2  in  Fig.  5.6).  The 
long  exposure  psf  for  the  screens  compensated  by  the  actual  tilt  is  shown  in  Fig.  5.11  (case 
1  in  Fig.  5.6).  The  long  exposure  psf  for  the  screens  compensated  by  the  phase/amplitude 
image  centroid  estimated  tilt  is  shown  in  Fig.  5.12  (case  3  in  Fig.  5.6).  Notice  how  the 
uncorrected  long  exposure  psf  is  shorter  and  wider  than  that  of  the  corrected  screens. 
Also,  notice  how  the  screen  compensated  by  using  the  phase/amplitude  centroid  estimate 
is  higher  and  narrower  than  the  screens  compensated  by  the  actual  average  tilt  and  the 
phase  centroid  estimate. 
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Figure  5.7  Point  spread  function  of  the  diffraction  limited  unaberated  screen. 
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Figure  5.8  Point  spread  function  (psf)  of  the  uncorrected  random  screen  for  the  short 
exposure  image  for  the  following  parameters  L/Lq  =  50,000,  D/vo  =■  10, 
Lol'f'o  =  200,  and  X/Lq  =  5.0  x  10“®. 


Figure  5.9  Point  spread  function  (psf)  of  the  uncorrected  random  screen  for  the  long 
exposure  image  for  the  following  parameters  L/Lo  =  50,000,  D/to  =  10, 
LoIto  =  200,  and  X/Lq  =  5.0  x  10“®. 
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Figure  5.10  Point  spread  function  (psf)  of  the  corrected  random  screen  using  the  image 
centroid  of  the  phase  screen  only  to  determine  the  tilt  for  the  long  exposure 
image  for  the  following  parameters  L/Lq  =  50, 000,  D/to  =  10,  Lo/to  =  200, 
and  A/Lo  =  5.0  x  10“®. 
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Figure  5.11  Point  spread  function  (psf)  of  the  corrected  random  screen  using  the  actual 
average  tilt  for  the  long  exposure  image  for  the  following  parameters  L/Lo^ 
50,000,  D/to  =  10,  Lo/ro  =  200,  and  X/Lo  =  5.0  x  10“®. 


Figure  5.12  Point  spread  function  (psf)  of  the  corrected  random  screen  using  the  image 
centroid  of  the  phase/amplitude  screen  for  the  long  exposure  image  for  the 
following  parameters  LfLo  =  50,000,  D/ro  =  10,  Lo/’’o  =  200,  and  A/Lq  = 
5.0  X  10"®. 


In  Fig.  5.13  the  strength  of  turbulence  is  held  constant  regardless  of  the  length  of 
propagation  allowing  for  an  isolation  of  the  effects  of  changing  LjLo.  For  longer  propa¬ 
gation  paths  the  relative  phase  perturbations  decrease  while  the  amplitude  perturbations 
increase.  Since  phase  perturbations  have  the  greatest  effect  on  the  psf,  the  Strehl  ratios 
rise  as  the  phase  perturbations  decrease. 

From  Fig.  5.13  one  can  see  using  the  amplitude/phase  screen  image  centroid  to 
estimate  tilt  correction  allows  for  a  higher  Strehl  ratio,  SRutaPj  over  the  Strehl  ratios, 
SRut  and  SRutp,  that  use  the  phase  to  estimate  the  correction.  Equations  (5.12)  and 
(5.13)  map  the  centroid  location  to  the  tilt  correction  required  to  move  the  centroid  to 
the  center  of  the  image.  Since  the  phase  as  well  as  the  amplitude  perturbations  affect  the 
position  of  the  centroid,  subtracting  the  average  phase  tilt  over  the  screen  will  not  take  into 
account  the  effect  the  amplitude  perturbations  have  on  the  image  centroid.  If,  however, 
the  phase  and  amplitude  screen  image  centroid  is  used  to  determine  the  tilt  correction,  the 
phase  and  amplitude  image  centroid  will  move  closer  to  the  center  of  the  image.  The  less 
movement  of  the  image  centroid  the  higher  the  Strehl  ratio  will  be  for  the  ensemble. 


Figure  5.13  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  LjLo  for  Djro  = 

10,  LqIto  =  200,  and  XjLo  =  5.0  x  10“®.  Strehl  ratios  are  plotted  for 
the  compensated  screen  with  phase/amplitude  centroid,  SRutap,  compen¬ 
sated  screen  with  actual  phase  tilt,  SRut,  compensated  screen  with  phase 
centroid,  SRutp,  and  uncompensated  screen,  SR. 
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In  Fig.  5.14,  as  the  strength  of  turbulence,  X-o/^oi  increases  the  uncompensated  Strehl 
ratio,  SR,  remains  relatively  constant.  As  Lo/to  increases  and  D/to  is  held  constant  the 
lower  order  phase  perturbations  increase  (i.e.,  more  tilt)  while  the  higher  order  pertur¬ 
bations  decrease  (i.e.,  flatter  screens).  The  mean  square  tilt  over  the  aperture  increases 
signiflcantly  as  seen  in  Fig.  5.3,  but  the  decrease  in  higher  order  phase  perturbations  coun¬ 
teracts  the  effects  of  the  larger  tilt.  Therefore,  as  the  screens  flatten  due  to  increases  in 
Lo/to,  the  Strehl  ratios  of  the  compensated  screens  will  be  higher.  As  before,  the  screens 
compensated  by  the  phase/amplitude  centroid  calculated  tilt  had  the  highest  Strehl  ratio. 


Figure  5.14  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  Lo/to  for  D/to  = 
10,  L/Lo  =  10000,  and  X/Lg  =  5.0  x  10~®.  Strehl  ratios  are  plotted  for  the 
compensated  screen  with  phase/amplitude  centroid,  SRutaP)  compensated 
screen  with  actual  phase  tilt,  SRut,  compensated  screen  with  phase  cen¬ 
troid,  SRutp,  and  uncompensated  screen,  SR. 
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Prom  Fig.  5.15  notice  how  increasing  A/Lo  has  nearly  the  same  effect  as  increasing 
L/Lq-  This  similar  relationship  is  due  to  the  mean  square  value  of  the  amplitude  and 
phase  perturbations’  dependence  upon  ^  as  seen  in  Figs.  2.5  and  2.6. 


Figure  5.15  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  A/Lo  for  Djro  = 
10,  L/Lq  =  10000,  and  Lo/to  =  200.  Strehl  ratios  are  plotted  for  the 
compensated  screen  with  phase/amplitude  centroid,  SRutaPi  compensated 
screen  with  actual  phase  tilt,  SRut,  compensated  screen  with  phase  cen¬ 
troid,  SRutp,  and  uncompensated  screen,  SR. 
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Figure  5.16  shows  how  increasing  the  ratio  D/tq  significantly  decreases  the  Strehl 
ratio.  Pried  parameter,  Tq,  can  be  interpreted  as  the  limiting  aperture  size  beyond  which 
the  turbulence  does  not  affect  the  resolution.  Therefore,  as  D/rg  gets  larger  the  more  the 
turbulence  limits  the  resolution.  The  improvement  in  the  Strehl  ratio  due  to  tilt  removal 
becomes  negligible  as  Djro  approaches  100.  Removing  the  tilt  from  the  screen  becomes 
negligible  since  as  the  aperture  diameter  increases  the  average  tilt  over  the  aperture  de¬ 
creases  significantly  as  shown  in  Fig.  5.5. 


Figure  5.16  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  D/to  for  X/Lg  = 
5.0  X  10“®,  LjLg  =  10000,  and  Lg/vg  =  200.  Strehl  ratios  are  plotted  for 
the  compensated  screen  with  phase/amplitude  centroid,  SRutaPi  compen¬ 
sated  screen  with  actual  phase  tilt,  SRut,  compensated  screen  with  phase 
centroid,  SRjjtPi  and  uncompensated  screen,  SR. 
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5.3  Comparisons  of  simulations  using  correlated  versus  uncorrelated  phase  and  am¬ 
plitude  screens 

In  this  section  the  same  calculations  are  implemented  as  in  the  last  section  but  this 
time  the  uncorrelated  phase  and  amplitude  screens  are  compared  to  calculations  using  the 
correlated  screens.  The  amplitude  and  phase  screens  for  the  uncorrelated  and  correlated 
case  will  both  be  calculated  using  Dp  =  1.5.  These  calculations  will  be  performed  to  deter¬ 
mine  if  implementing  the  cross  correlation  affects  the  results  in  any  way.  The  parameters, 
L/Lo  and  Lo/ro,  were  chosen  to  be  studied  since  varying  these  parameters  had  different 
effects  on  the  Tmse  and  the  Strehl  ratios. 

5.3.1  Mean  Square  Error  in  the  Tilt  Estimation  with  correlated  versus  uncorre¬ 
lated  screens.  The  mean  square  error  in  the  tilt  estimation  was  calculated  for  both  the 
correlated  and  uncorrelated  phase  and  amplitude  screens  to  determine  whether  the  cross 
correlation  has  any  effect  on  the  results.  Figures  5.17  and  5.18  show  how  Tms  and  Tmae 
vary  with  LjLo  for  the  correlated  and  uncorrelated  screens,  respectively.  The  uncorre¬ 
lated  case  seemed  to  have  a  slightly  higher  Tms  and  Tmse  for  the  amplitude/phase  screen. 
Figures  5.19  and  5.20  show  how  Tms  and  Tmse  vary  with  Lofro  for  the  correlated  and 
uncorrelated  screens,  respectively.  This  time,  the  correlated  screens  case  seemed  to  have 
a  slightly  higher  Tmsi  but  slightly  lower  Tmse  for  the  amplitude/phase  screen.  As  one  can 
see  the  mean  square  error  in  the  tilt  estimation  is  not  significantly  affected  by  the  cross 
correlation  of  the  amplitude  and  phase. 

5.3.2  Strehl  Ratio  Calculation  with  correlated  versus  uncorrelated  screens. 
Strehl  ratios  were  calculated  for  both  the  correlated  and  uncorrelated  phase  and  amplitude 
screens  to  determine  whether  the  cross  correlation  has  any  effect  on  the  results.  Figures 
5.21  and  5.22  show  how  the  Strehl  ratios  vary  with  LjLo  for  the  correlated  and  uncorrelated 
screens,  respectively.  Figures  5.23  and  5.24  show  how  the  Strehl  ratios  vary  with  Lofro 
for  correlated  and  uncorrelated  screens,  respectively.  The  Strehl  ratios  for  the  correlated 
and  uncorrelated  screens  are  just  slightly  different,  but  SRut  and  SRjjtp  are  consistently 
slightly  higher  for  the  correlated  screens  than  the  uncorrelated  screens. 
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Figure  5.17 


Figure  5.18 


The  amplitude  and  phase  screen  are  correlated  by  the  cross  correlation. 
The  mean  square  error  in  the  tilt  estimation,  T^nsej  and  the  mean  square 
tilt,  Tms,  are  plotted  versus  LfLo  for  Dfro  =  10,  Z-oAo  =  200,  and  A/Lo  = 
5.0  X  10“®.  The  is  calculated  twice,  once  using  the  phase  and  amplitude 
screen  case  and  once  for  the  phase  screen  only  case.  The  Tms  is  calculated 
by  computing  the  average  slope  of  the  phase  over  the  screen.  Dp/Lo  = 
D^/Lo  =  DpLo  =  D^/Lo  =  1.5. 


The  amplitude  and  phase  screen  are  uncorrelated.  The  mean  square  error 
in  the  tilt  estimation,  Tmaa  and  the  mean  square  tilt,  Tms,  are  plotted  versus 
LfLo  for  D/to  =  10,  Lq/to  =  200,  and  A/Lp  =  5.0  x  10~®.  The  Tmse  is 
calculated  twice,  once  using  the  phase  and  amplitude  screen  case  and  once 
for  the  phase  screen  only  case.  The  Tms  is  calculated  by  computing  the 
average  slope  of  the  phase  over  the  screen.  D^ILq  =  1.5  and  D^/Lo  =  1.5. 
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Figure  5.19 


Figure  5.20 


ro 


The  amplitude  and  phase  screen  are  correlated  by  the  cross  correlation. 
The  mean  square  error  in  the  tilt  estimation,  Tmsef  and  the  mean  square 
tilt,  Tffis,  are  plotted  versus  Lofro  for  Djro  =  10,  L/Lq  =  10000,  and 
\/Lo  =  5.0  X  10“®.  The  T,nse  is  calculated  twice,  once  using  the  phase  and 
amplitude  screen  case  and  once  for  the  phase  screen  only  case.  The  Tms 
is  calculated  by  computing  the  average  slope  of  the  phase  over  the  screen. 
DpiLo  =  D^ILo  =  DpLo  =  D^lLo  =  1.5. 


The  amplitude  and  phase  screen  are  uncorrelated.  The  mean  square  error 
in  the  tilt  estimation,  Tmse^  and  the  mean  square  tilt,  Tms,  are  plotted  versus 
LoIto  for  D/to  =  10,  LjLo  =  10000,  and  A/Lo  =  5.0  x  10“®.  The  Tmse  is 
calculated  twice,  once  using  the  phase  and  amplitude  screen  case  and  once 
for  the  phase  screen  only  case.  The  Tms  is  calculated  by  computing  the 
average  slope  of  the  phase  over  the  screen.  D^/Lq  =  1.5  and  D^fLg  =  1.5. 
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Figure  5.21 


Figure  5.22 


The  amplitude  and  phase  screen  are  correlated  by  the  cross  correla¬ 
tion.  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  LfLo 
for  D/to  =  10,  Lo/tq  =  200,  and  X/Lo  =  5.0  x  10“®.  Strehl  ra¬ 
tios  are  plotted  for  the  compensated  screen  with  phase/amplitude  cen¬ 
troid,  SRutaPi  compensated  screen  with  actual  phase  tilt,  SRut,  com¬ 
pensated  screen  with  phase  centroid,  SRutp,  and  uncompensated  screen, 
SR.  Dp/Lo  =  D^/Lo  =  DpLo  =  D^lLo  =  1.5. 


The  amplitude  and  phase  screen  are  uncorrelated.  Strehl  ratios  of  the 
simulated  wavefront  are  plotted  versus  L/Lo  for  D/to  =  10,  Lo/vo  =  200, 
and  X/Lo  =  5.0  x  10“®.  Strehl  ratios  are  plotted  for  the  compensated  screen 
with  phase/amplitude  centroid,  SRutaPi  compensated  screen  with  actual 
phase  tilt,  SRut,  compensated  screen  with  phase  centroid,  SRutp,  and 
uncompensated  screen,  SR.  D^fLo  =  1.5  and  D^fLo  =  1.5. 


5-25 


Figure  5.23 


Figure  5.24 


The  amplitude  and  phase  screen  are  correlated  by  the  cross  correla¬ 
tion.  Strehl  ratios  of  the  simulated  wavefront  are  plotted  versus  LoItq 
for  Djro  =  10,  LjLo  =  10000,  and  A/Lq  =  5.0  x  10~®.  Strehl  ra¬ 
tios  are  plotted  for  the  compensated  screen  with  phase/amplitude  cen¬ 
troid,  SRutap,  compensated  screen  with  actual  phase  tilt,  SRut,  com¬ 
pensated  screen  with  phase  centroid,  SRutp^  and  uncompensated  screen, 
SR.  Dp/Lo  =  Dj/Lo  =  Df/Lo  =  D^lLo  =  1.5. 


The  amplitude  and  phase  screen  are  uncorrelated.  Strehl  ratios  of  the 
simulated  wavefront  are  plotted  versus  Lojro  for  D Jr „  =  10,  L/Lq  =  10000, 
and  A/Lo  =  5.0  x  10~®.  Strehl  ratios  are  plotted  for  the  compensated  screen 
with  phase/amplitude  centroid,  SRutap>  compensated  screen  with  actual 
phase  tilt,  SRut,  compensated  screen  with  phase  centroid,  SRutPj  and 
uncompensated  screen,  SR.  D^/Lq  =  1.5  and  D^/Lq  =  1.5. 
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5.3.3  Conclusions  about  using  correlated  versus  uncorrelated  screens.  In  this 
study  uncorrelated  screens  were  used  since  the  amplitude  perturbations  were  not  suffi¬ 
ciently  represented  in  the  FS  expansion.  There  are  slight  differences  between  using  cor¬ 
related  versus  uncorrelated  screens,  but  the  reduction  in  the  amplitude  perturbations  due 
to  the  under  representation  proved  to  be  of  greater  concern.  Correlated  screens  could 
have  been  used  in  this  study  if  more  FS  coefficients  were  used  in  the  amplitude  screen 
generation,  but  the  matrix  dimensions  required  (1080  x  1080)  made  the  calculation  too 
computationally  intensive. 
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VI.  Conclusions  and  Recommendations 

6.1  Conclusions 

A  new  atmospheric  turbulence  screen  generator  was  developed  for  use  in  performance 
simulations  of  adaptive  optics  algorithms  valid  over  a  wide  range  of  atmospheric  turbulence 
parameters,  incorporating  both  phase  and  amplitude  effects.  Most  screen  generators  only 
calculate  the  phase  perturbations  over  the  aperture  [2,4,11,12,15].  The  screen  generator 
developed  here  accounts  for  diffraction  effects  caused  by  weak  turbulence  by  incorporating 
the  phase,  amplitude,  and  cross  statistics  of  the  perturbed  optical  field.  The  wavefront’s 
phase  and  amplitude  pertmbations  were  taken  from  the  correlation  functions  developed 
by  Lee  and  Harp.  The  cross  statistics  of  Lee  and  Harp’s  weak  turbulence  model  were 
developed  here  for  the  first  time  and  were  incorporated  into  the  screen  generator.  Since 
the  amplitude  perturbations  cannot  be  corrected  with  adaptive  optics,  the  theoretical 
limit  of  an  adaptive  optics  system  is  bounded  by  the  effects  of  the  amplitude  pertubations. 
Using  both  the  amplitude  and  phase  perturbations,  one  can  more  accurately  simulate  the 
performance  of  an  adaptive  optics  system.  Since  a  weak  index  of  refraction  fluctuation 
model  was  used  in  the  development,  the  screen  generator  will  not  predict  effects  of  strong 
turbulence,  such  as  the  satmation  of  the  log-amplitude  fluctuations. 

For  independent  uncorrelated  phase  and  amplitude  screens  the  mean  square  value  and 
structure  functions  of  the  Monte  Carlo  simulation  were  within  1%  of  the  theoretical  values. 
Due  to  difierences  in  the  amplitude  and  pheise  power  spectral  densities  and  computational 
limits  on  the  matrix  size  in  the  calculations,  the  simulations  for  the  correlated  screens  were 
within  7%  of  the  theoretical  values  for  the  amplitude,  1%  for  the  phase,  and  1%  for  the 
cross  statistics. 

Incorporating  the  cross  correlation  into  the  simulations  had  a  negligible  effect  on 
the  results  obtained  and  therefore  the  uncorrelated  phase  and  amplitude  screens  were 
used  to  more  accurately  represent  the  phase  and  amplitude  perturbations.  Monte  Carlo 
experiments  performed  using  the  screen  generator  showed  the  amplitude  perturbations 
can  significantly  reduce  the  accuracy  of  full-apertmre  tilt  estimation  using  image  centroid 
motion. 
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In  contrast,  using  the  same  amplitude/phase  screen  image  centroid  to  estimate  tilt 
correction  allows  for  a  higher  Strehl  ratio  than  using  only  the  phase  to  estimate  the  correc¬ 
tion.  The  tilt  estimation  algorithm  maps  the  image  centroid  location  to  the  tilt  correction 
required  to  move  the  centroid  to  the  center  of  the  image.  Since  the  phase  as  well  as  the 
amplitude  perturbations  affect  the  position  of  the  centroid,  subtracting  the  average  phase 
tilt  over  the  screen  will  not  take  into  account  the  effect  the  amplitude  pertmrbations  have 
on  the  image  centroid.  If,  however,  the  phase  and  amplitude  screen  image  centroid  is 
used  to  determine  the  tilt  correction,  the  phase  and  amplitude  image  centroid  will  move 
closer  to  the  center  of  the  image.  The  less  movement  of  the  image  centroid  the  higher  the 
Strehl  ratio  will  be  for  the  long  exposure  image.  Therefore,  using  the  tilt  estimator  on  the 
combined  phase  and  amplitude  screen  resulted  in  higher  Strehl  ratios  over  subtracting  the 
average  tilt  over  the  aperture. 

6.2  Recommendations  for  Further  Study 

6.2.1  Screen  Generation.  In  this  thesis  the  turbulence  is  modeled  as  uniform 
over  the  propagation  path  which  is  often  appropriate  for  horizontal  propagation  paths.  A 
layered  approach  could  be  implemented  to  accurately  represent  the  strength  of  turbulence 
profile  up  through  the  atmosphere.  Using  Taylor’s  frozen  flow  hypothesis  one  could  also 
incorporate  the  time  statistics  into  the  amplitude  and  phase  screen  generator.  As  in 
Welsh’s  paper  on  phase  screen  generation  one  could  incorporate  the  layered  representation 
of  the  phase  and  amplitude,  incorporate  the  temporal  statistics,  and  provide  for  off-axis 
calculations  to  include  the  effects  of  anisoplanatism  [15]. 

6.2.2  Simulations.  The  atmospheric  turbulence  generator  developed  in  this  thesis 
can  be  used  to  study  the  effects  of  amplitude  perturbations  on  all  sorts  of  adaptive  optics 
alogorithms.  The  screen  generator  can  also  be  used  to  determine  the  effects  of  different 
atmospheric  conditions  on  imaging  systems  or  laser  propagation.  Including  the  phase  and 
amplitude  pertmrbations  in  simulations  will  allow  one  to  more  accurately  simulate  the 
complete  performance  of  an  adaptive  optics  system. 
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Monte  Carlo  experiments  performed  using  the  screen  generator  showed  the  amplitude  perturbations  can  significantly  reduce 
the  accuracy  of  full-aperture  tilt  estimation  using  image  centroid  motion.  However,  since  the  amplitude  perturbations  affect 
the  image  centroid,  the  tilt  estimate  using  the  combined  amplitude  and  phase  screen  allows  for  a  higher  Strehl  ratio  than 
using  only  the  phase  to  estimate  the  correction. 
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